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Abstract 



In this thesis we examine the method of counting BB events produced in the 
BaBAR experiment. The original method was proposed in 2000, but improve- 
ments to track reconstruction and our understanding of the detector since that 
date make it appropriate to revisit the B Counting method. We propose a new 
set of cuts designed to minimize the sensitivity to time-varying backgrounds. 
We find the new method counts BB events with an associated systematic un- 
certainty of ±0.6%. 
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Chapter 1 



Introduction 



The BaBAR experiment is located at the Stanford Linear Accelerator Center 
(SLAC). Data from e+e" annihilations at a centre-of-mass (CM) energy ^/s « 
10.58 GeV are taken at the SLAC PEP-II storage-rings with the BaBAR detector. 
This chapter provides a brief introduction to the experiment. 

PEP-II is a 'i?-Factory' — it is designed to produce a high number of B 
mesons. This is achieved through the process e+e~ T{AS) — > BB. Beams of 
9.0 GeV electrons and 3.1 GeV positrons collide with a CM energy of 10.58 GeV, 
which is the peak of the T(4S') resonance. Running at this energy is desirable 
because the T(4S') meson (a bound bb state) decays more than 96% of the time 
to a BB state. The design luminosity of PEP-II (3 x 10'^^ cm^s~^) was exceeded 
in 2001 and by 2006, peak luminosities above 12 x 10'^^ cm^s~^ were recorded. 



PEP II 

Low Energy 




•m 3 km N- 

Figure 1.1: Diagram of PEP-II linac and storage-rings. 

Asymmetric beam energies cause the decay products to have a Lorentz boost 
of P"f = 0.56 relative to the laboratory. This is one of the most important design 
features of the BaBAR experiment. 

In the rest frame of each T(45) particle, the B mesons are created almost at 
rest. By boosting along the z axis, the B mesons travel a measurable distance 
inside the detector before decaying. The mean lifetime of a -B meson is very 
short (around 1.5 ps), but by reconstructing the decay vertices of each B it 
possible to calculate how far each travelled and hence the lifetimes. 

For some analyses (such as studies of CP violation in the BB system) the 
relative difference in decay times is especially important. The BaBAR detector 
is specifically designed to enable precise measurements of this quantity. 
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PEP-II is run at the T(4S') resonance most of the time, but it is important 
to occasionally run 'off-peak'. Typically just over 10% of the time, PEP-II is 
run around 40 MeV below the resonance peak. This is achieved by lowering the 
energy of the electron beam, which reduces the boost by less than one percent. 
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Figure 1.2: The spectrum of hadron production near PEP-II's operational CM 
energy. The curve shows the cross-section for inclusive production 
of hadrons (vertical axis) as a function of CM energy. The peak at 
the T(4S') resonance is clearly visible. The plot is originally from 
Upsilon Spectroscopy by Besson and Skwarnicki [1]. 



Continuum is the name used to describe all non-BB events produced in 
the detector. The continuum contains many types of events, by far the most 
common of which is e^e^ — )■ e^e^ (Bhabha events) which have a cross-section 
close to 40 nb. For each type of quark q (apart from t), events of the type 
e^e~ — > qq have a cross-section of order 1 nb. The other leptons (/i+^^ or 
r+T^) are produced with cross-sections also close to 1 nb. 

These cross-sections all scale in a known way when PEP-II runs at the de- 
creased off-peak energy, while T{AS) (and hence BB) production is 'switched 
off'. The off-peak data can then be used to understand on-peak backgrounds 
and when scaled in the correct way can be used to count B events in on-peak 
data. Any B physics analysis at BaBAR will make use of the off-peak data, and 
it is an integral part of B Counting. 
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1.1 CP Violation 



The BaBAR experiment was designed to achieve a number of physics goals. The 
primary motivation for the experiment was to study CP violation in the decays 
of B mesons [2]. A significant result '3] demonstrating the existence of CP 
violation in this sector was published in 2001, two years after running began. 

CP is the product of two quantum mechanical operators. Charge conjuga- 
tion (C) interchanges particles with their anti-particles and parity (P) changes 
the 'handedness' of a co-ordinate system, i.e. P sends (i,x) (<, — x). 

Independently, C and P are not symmetries of nature and in particular, 
neither is a symmetry of the weak nuclear force. For example, left-handed 
neutrinos are transformed to right-handed neutrinos by the operation of P, 
and to left-handed anti-neutrinos by C. Neither of these states is observed to 
participate in the weak interaction. 

Cronin and Fitch were awarded the 1980 Nobel Prize in Physics for demon- 
strating the combination CP was also not a symmetry of nature They 
observed this in the decay of neutral kaons, specifically through the mode 
K'l — > TT+vr^ which would not occur if CP were conserved. Following this, 
CP violation was also anticipated in B meson decays, and was expected to be 
larger and more varied than in the kaon system. Investigating this was one of 
the main reasons for running the BaBAR experiment. 

There are several reasons for studying CP violation with an experiment 
like BaBAR. Measurements of CP violation arc important to test the Standard 
Model (SM) of particle physics. Any results outside SM predictions by defi- 
nition indicate new physics. A second reason relates to the matter- antimatter 
asymmetry of the universe. CP violation is one of the conditions proposed by 
Sakharov [5] to explain how the observed asymmetry could occur. 

In the Standard Model, CP violation is governed by the weak quark mixing 
matrix, known as the Cabibbo-Kobayashi-Maskawa (CKM) matrix [51 [7] • This 
takes the form: 



where each Vij is complex, and the probability of a transition between two types 
of quark i and j is proportional to p. Here u, d, s, c, t and b represent up, 
down, strange, charm, top and bottom quarks respectively. 

By studying a variety of different decay channels, BaBAR can be used to 
measure or put limits on several components of the CKM matrix. In particular, 
many channels seen at BaBAR allow the parameter 



to be measured. Other searches for CP-violation, such as in gluonic penguin de- 
cays allow the measurement of different CKM matrix elements (see for example 




(1.1) 




(1.2) 
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[S]) including: 



7 = Arg 



( 




) 



(1.3) 



1.2 Other Physics at BaBAR 



The BaBAR experiment has muhiple uses apart from searching for CP violation 
in B decays. It is a high luminosity experiment, and as well as producing many 
B mesons (enabling study of rare B processes), collisions create many other 
species including charmonium states (such as J/ip ) , other charmed mesons (such 
as D), T leptons and two- photon states. 

Data from BaBAR are used to study each of these to refine and test the 
Standard Model. The cross-section for cc production in the BABAR detector is 
slightly larger than for bb, while the r+r" cross-section is slightly smaller. 

The charm physics research at BaBAR includes measuring D meson and 
charm baryon lifetimes, rare D meson decay and D mixing. An important 
result demonstrating evidence for D'^D^ mixing observed in Z?" Kn decays 
was published in 2007 [9J. 

An analysis group is specifically devoted to studying charmonium (cc) states. 
This group studies several types of physics. This includes initial state radiation 
(ISR) production (such as e~^e~' — > 7/5fl7r+7r~ J/V'), fully reconstructed B events 
with charmonium final states (for example B J/ipX or B tp{2S)X where X 
represents other unspecified particles) and studies of exotic states like Y(4260). 

Other groups study charmless hadronic decays of B mesons (e.g. B tttt), 
leptonic c and b decays (e.g. D,B ii^i, where £ = e, fj, or t) and radiative 
penguins. The last of these describes processes involving radiative loops such 
as 6 — *■ 57. Radiative penguin events provide one way to test theories beyond 
the Standard Model because new virtual particles could appear in the loop with 
detectable effects. 

Also, research at BABAR includes measurements of the properties of the r 
lepton produced through e+e~ t^t^. An important area of research in 
the r sector is the search for lepton flavour violation through decays such as 

— » or T+ — » e~^e~e~^ . 
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The BaBAR Detector 



This chapter provides a brief introduction to the hardware of the BaBAR de- 
tector. Detailed descriptions of the detector [1^ and also of PEP-II [13] are 
available elsewhere. 

It will be useful to briefly explain the (right-handed) co-ordinate system used 
in the BaBAR experiment. 

• The ~\-z axis is the direction of the high-energy (electron) beam. 

• The +y axis is vertically upwards. 

• The azimuthal angle (fi lies in the xy plane. It is zero on the +x axis, and 
increases towards the +y axis. 

• The polar angle 9 is measured from the z axis. It is zero on the +z axis 
and TT on the —z axis. 

The beam energy asymmetry of PEP-II means the optimal detector shape is 
also asymmetric. A schematic of the detector is shown in Fig. 12.11 Full details 
of the detector and its sub-systems can be found in ^14] and [12] . The detector 
was designed to meet a large number of requirements which include: 

• performing tracking on particles with transverse momentum (pt) between 
- 60 MeV/c and ~ 4 GeV/c; 

• detecting photons and neutral pions with energies between ^ 20 MeV and 
~ 5 GeV, and; 

• distinguishing electrons, muons, pions, kaons and protons over a wide 
range of momenta. 

To meet these requirements, the detector was constructed with five major 
subsystems. The innermost components are a five-layer double-sided silicon 
vertex tracker (SVT), which provides precision tracking information and a 40- 
layer drift chamber (DCH), which provides momentum information for charged 
tracks. Outside these is a particle identification system, a detector of internally- 
reflected Cerenkov light (DIRC). Surrounding all these is an electromagnetic 
calorimeter (EMC) built from more than six thousand CsI(Tl) crystals. The 
EMC records energy deposited by electrically charged particles and photons 
and also has some sensitivity to neutral hadrons such as neutrons and K^. The 
EMC is able to detect neutral pions by recording the photons from tt*^ 77. 
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Finally, outside the superconducting coil is the muon detection system, called 
the Instrumented Flux Return (IFR) consisting mostly of Limited Streamer 
Tube (LST) modules. 




3-2001 
8553 A50 



Figure 2.1: Cross-section of the BaBar detector. 



2.1 The Silicon Vertex Tracker 

The silicon vertex tracker (SVT) is the component of the BaBAR detector closest 
to the interaction point (IP). It provides the most precise angular measurements 
of tracks because further from the IP, track precision is limited by multiple 
scattering. The SVT consists of five cylindrical layers of double-sided silicon- 
strip detectors. 

In total, there are 340 silicon detectors within the SVT, with around 150,000 
readout channels. The individual silicon detectors vary in size between 43 x 42 
mm^ (z X (p) and 43 x 42 mm^. Each detector is double-sided, with strips to 
measure (j) running parallel to the beam axis (</) strips) and strips to measure z 
{z strips) running orthogonally on the opposite side. The total area covered by 
the five layers of silicon detectors is around one square metre. 

The ends of each strip are held at a potential difference of around 35-45 V. 
When a charged particle passes through the silicon, it creates free electrons and 
electron-hole pairs. This charge is collected by the applied voltage, and after 
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being amplified can be read out as a 'hit' signal. The track position is calculated 
by combining all the hits from individual silicon detectors. 

The resolution of the SVT depends on track angle and varies between layers. 
For tracks at normal incidence, the resolution of inner layers is approximately 
15 ^m in z and 10 /xm in </>. For outer layers the resolutions are approximately 
35 /Ltm and 20 /zm respectively. 

Some very low momentum particles such as slow pions in D* decay have in- 
sufficient energy to reach outer components of the detector. For these particles, 
the SVT provides the only tracking information. 

2.2 The Drift Chamber 

The multi-wire drift chamber (DCH) provides precise tracking information and 
to a lesser extent is used to perform some particle identification (PID). It was 
constructed on the campus of the University of British Columbia at TRIUMF. 
The DCH is cylindrical and surrounds the SVT. It is filled with a gas mixture 
of helium and isobutane, which is kept 4 mb above atmospheric pressure. 

The DCH is strung with nearly thirty thousand wires arranged into 40 layers 
(and ten 'super layers') of hexagonal cells. The layout of these cells for 16 layers 
is shown in Fig. 12.21 A tungsten-rhenium 'sense' wire in the centre of each cell 
is held at a positive high voltage, of around 1960 V. The surrounding gold- 
coated aluminium 'field' wires are at ground potential. Cells on the boundary 
of a superlayer have two gold-coated aluminium 'guard' wires. These are held at 
340 V and improve the performance of the cells and ensure the gain of boundary 
cells is the same as that of inner ones. 

Of the ten superlayers, four are 'stereo'. The wires of these layers are strung 
at a slight angle (between ±45 and ±76 mrad) relative to the z axis. This is 
important for tracking as it enables 9 co-ordinates to be measured. 

A charged particle moving through the drift chamber ionizes the gas. The 
ionized gas molecules within a cell move towards a ground wire and the electrons 
released from these molecules are accelerated towards the sense wire. In turn 
these electrons collide with and ionize other gas molecules. This causes an 
'avalanche' of negative charge (up to a gain of 50,000) to arrive at the sense 
wire. For a particle energetic enough to exit the drift chamber, up to 40 DCH 
'hits' of this sort (one per cell) can be recorded. 

The drift chamber is contained within the 1.5 T field of the super-conducting 
coil. Measurements of track curvature in the drift chamber allow momenta to 
be determined. The momentum resolution of the drift chamber for a particle of 
transverse momentum pt (i.e. the component of momentum perpendicular to 
the z axis) in units of GeV/c is: 

ap, = (0.13 ± 0.01)% • Pt + (0.45 ± 0.03)%. (2.1) 

To be measured in the drift chamber (i.e. to travel past the inner radius of the 
DCH and trigger enough wires to provide a good measurement), a track must 
have transverse momenta of at least 100 MeV/c. 
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Figure 2.2: A schematic drawing of 16 layers of DCH drift cells. The lines con- 
necting wires are included for visualisation purposes only. Numbers 
in the 'Stereo' column indicate the angle in mrad of the stereo layers. 



Different species of particles have a characteristic rate of energy loss with 
distance (dE/dx), and the drift chamber can provide PID for low-momentum 
tracks by measuring this as a particle moves through the chamber. The dE/dx 
distributions for different particles as a function of momentum are shown in 
Fig. 12.31 Particle ID is achieved by recording the total charge deposited in each 
cell. Corrections are made to account for several sources of bias including signal 
saturation, variations in gas pressure and differences in cell geometry. 

For low-momentum tracks (i.e. those below 1 GeV/c), the overall dE/dx 
resolution is around 7%. The final PID combines this information with that 
gathered from the DIRC and other subsystems. 
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Figure 2.3: Drift Chamber dE/dx as a function of particle momentum. The 
scatter points are from beam scan data and the curves are appro- 
priately parameterized Bethe-Bloch curves. The letters have their 
usual meanings, labelling the curves muon, electron, pion, kaon, 
proton and deuteron respectively. 



2.3 The DIRC 

The function of the detector of internally-reflected Cerenkov light (DIRC) is to 
perform particle identification. In particular, it especially separates kaons and 
pions of momenta between 2.0 and 4.0 GeV/c. It records the Cerenkov radiation 
of charged particles passing through a ring of fused-silica bars located outside 
the drift chamber. This type of light is emitted whenever a charged particle 
passes through a medium with a velocity greater than the speed of light in that 
medium. 

For a medium of refractive index n (for fused-silica, n « 1.47), a particle 
moving with velocity v = (3c (where c denotes the speed of light in a vacuum) 
emits a characteristic cone of Cerenkov light at an angle cos(6'c) — l/(/3n). 
When the Cerenkov angle 9c of a particle is measured in the DIRC, momentum 
information from the drift chamber allows the mass (and hence type) of the 
particle to be determined. 

Ring-imaging Cerenkov counters have commonly been used in particle de- 
tectors since the 1980s, but the BaBAR DIRC differs from all previous ones by 
utilising total internal reflection. The DIRC radiator comprises 144 synthetic 
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Figure 2.4: A schematic drawing of a DIRC radiator bar (left) and imaging 
region. 



quartz bars of length 4.9m which function both as light guides and as a medium 
for Cerenkov radiation. As shown in Fig. 12.41 light is carried towards DIRC 
electronics in the imaging region — located at the backward-end of the detec- 
tor — by successive internal reflections. A mirror is placed at the forward-end 
of the radiator bars to ensure forward-travelling photons are reflected back to 
arrive at the DIRC electronics. 

At the DIRC imaging region, the light expands into the 'Standoff Box' which 
is filled with six cubic metres of purified water of a similar refractive index to 
the quartz. The photons are detected at the rear of this by almost 11,000 
photomultiplier tubes (PMTs), each of diameter 2.82 cm. The position and 
arrival-time of PMT signals are then used to calculate the Cerenkov angle of 
particles passing through the DIRC radiator. 

At the high-momentum end of the DIRC's functional window (around 4.0 
GeV/c), the difference in the Cerenkov angle of pions and kaons is 3.5 mr. At 
this momentum, the DIRC is able to separate the two species to around 3cr. 
This increases for lower momenta, up to a separation of as much as lOcr at 2.0 
GeV/c. 

2.4 The Electromagnetic Calorimeter 

The Electromagnetic Calorimeter (EMC) measures electromagnetic showers in 
the energy range 20 MeV to 9 GeV. It comprises 6580 thallium-doped caesium 
iodide (CsI(Tl)) crystals which sit in 56 rings in a cylindrical barrel and forward 
end-cap which together cover a solid angle of —0.775 < cos6' < 0.962 in the 
laboratory frame. This coverage region was chosen because few photons travel 
in the extreme backward direction. The barrel part of the EMC is located 
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between the DCH and the magnet cryostat. 

Thallium- doped Csl is used because of its favourable properties. It has a 
high light yield, a small Moliere radius and a short radiation length. These 
allow excellent resolution of both energy and angle within the compact design 
of the detector. 



Figure 2.5: A cross-section of the EMC (top half). The EMC is symmetric in 



Altogether, the EMC weighs more than 26 tonnes. The crystals are held 
from the rear by an aluminium support system. Similarly, cooling equipment 
and cables are also located at the rear. This configuration minimises the amount 
of material between the interaction point (IP) and the crystals. The support 
structure is in turn attached to the coil cryostat. The layout of the aluminium 
support and the 56 crystal rings can be seen in Fig. 12.51 

High-energy electrons and photons create electromagnetic showers when 
travelling through the EMC crystals. The shower size depends on the prop- 
erties of the particle. The crystal absorbs the shower photons, and re-emits 
them as visible scintillation light. When tested with a source of 1.836 MeV 
photons, the average light yield per crystal is 7300 emitted photons per MeV of 
energy. The scintillation photons are read out with a package consisting of two 
silicon photodiodes and two preamplifiers mounted on each crystal. 

There is a degradation in time due to radiation damage inside the EMC, and 
the crystals are constantly monitored to gauge their individual performance. 
This is achieved by recording the EMC response both to known radioactive 
sources and to real events (i.e. from colliding beams). The latter involves cali- 
brating the system by measuring particles of a known type, energy and position 
with the EMC. A variety of events are used for this calibration, including e+e^ 
(Bhabhas), 77 and 




(f). Dimensions are given in mm. 
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2.5 The Muon Detection System 

The IFR consists of three sections — a central cyhndrical barrel and two end 
caps. When BaBAR was commissioned, the IFR was filled with Resistive Plate 
Counters (RPCs). Due to degradation and diminishing efficiency over the course 
of running, it became necessary to replace the barrel RPCs between 2004 and 
2006 with a different technology (LSTs). RPCs are still used in the end caps. 

The IFR is designed to detect deeply penetrating particles such as muons 
and neutral hadrons (mostly and neutrons). The system covers a large solid 
angle and was designed to have high efficiency and high background rejection. 

The IFR structure, shown in Fig. 12.61 is made from large quantities of iron, 
segmented into 18 plates of thickness varying between 2cm (nearest the beam) 
and 10cm (the outermost layer). The total thickness of the barrel is 65 cm, 
and the end caps 60 cm. The iron (which also functions as the superconducting 
magnet's flux return) is optimized to be a muon filter and neutral hadron ab- 
sorber allowing the best identification of muons and K^s. The iron layers are 
separated by a 3.2 cm gap where up to 19 RPCs, LSTs, or layers of brass can 
be housed. 



Figure 2.6: The IFR barrel and forward (FW) and backward (BW) endcaps. 

RPCs and LSTs are different designs which achieve the same goal — to 
provide accurately timed, precise two-dimensional positions of charges passing 
through. An RPC is filled with a gas mixture of argon, Freon and isobutane. 
Two graphite-coated surfaces located on either side of the gas are held at a 
potential difference of 6700-7600 V. Ionising particles create a streamer between 
these two surfaces and their presence is read-out by pairs of capacitive strips 
running parallel and perpendicular to the beam. 

Linseed oil was used to coat inner surfaces, and degradation of this and the 
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graphite surfaces contributed to efficiency losses as the experiment progressed. 
By 2002, the muon identification efficiency had dropped from 87% (at the start 
of running) to 78% at a pion misidentification rate of 4%. For this reason, 
LSTs were installed to replace all RPCs in the barrel. Two sextants were fully 
replaced in summer 2004, and the remaining four sextants in summer 2006. The 
RPCs in the forward end cap were replaced with new, more robust RPCs in fall 
2000 [H]. 

LSTs function slightly differently from RPCs. An LST consists of either 
seven or eight gas-filled cells. Through the length of each cell runs a wire held 
at high voltage (HV) of around 5500 V. To stop the wire sagging, several wire- 
holders are fitted along the cell's length. The gas mixture in the LST cells is 
predominantly carbon dioxide, and also contains smaller amounts of isobutane 
and argon. 

A charged particle passing through the cell ionizes the gas and creates a 
streamer which can be read out from the wire. This provides the (/> co-ordinate. 
The streamer also induces a charge on a plane below the wire. Running parallel 
to the wires are a series of conducting strips (z-planes) which detect this charge 
and provide the z co-ordinate. 

For the upgrades in 2004 and 2005, the nineteen layers of RPCs in each barrel 
sextant were replaced by twelve layers of LSTs and six of brass. The outermost 
layer of RPCs was inaccessible, so no LSTs could be installed there. The brass 
was installed in every second layer starting with the fifth both to compensate 
for the loss of absorption between the two outer layers and to increase the total 
absorption length. 
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Chapter 3 



Particle Lists 



The BaBAR database of events (including more than 400 mihion BE events) is 
very large and several stages of processing exist to reduce the time and storage 
requirements of analysts. The final selected data are stored in ntnples, and 
analysts must make a number of decisions regarding the types of information 
that should be stored. For example, it is often desirable to store only events of 
a certain type, and necessary to define what properties a track candidate should 
have before it is considered relevant. 

As well as containing well-measured particles (charged tracks and neutrals) 
from genuine physics events, raw data includes a lot of other ('junk') information 
which must be filtered as much as possible. This includes particles from cosmic 
rays, and beam-gas and beam-wall interactions. The software used by the BaBar 
Collaboration is flexible enough to give analysts control over this filtering. 

The term particle list means a list of all charged or neutral candidates in 
a particular event which pass certain criteria. These criteria are varied and 
depending on the type of list can include limits on track angle, momentum, 
distance of the track from the interaction point and many other quantities. A 
candidate is a collection of hits within the detector (for neutral particles this 
takes the form of energy deposits in certain detector components) which after 
reconstruction have been identified as a likely particle. Candidates can be either 
charged or neutral, and there is no guarantee a candidate is in fact anything 
other than background noise — it only has to appear to be a genuine particle. 

A charged candidate is a track associated with a particle moving inside the 
detector which (because of its charge) curves in the magnetic field. A neutral 
candidate is any EMC cluster which has no matching track. 

3.1 Lists Used in B Counting 

The most basic lists, on which others are based are ChargedTracks and CalorNeu- 
tral. These contain every charged and neutral candidate after event reconstruc- 
tion. 

The lists used in the past for B Counting are ChargedTracksAcc (for charged 
candidates) and GoodNeutralLooseAcc (neutrals). There is also an additional 
requirement based on the number of charged candidates on the GoodTracksAc- 
cLoose list. 

The lists used in this analysis are slightly different to these and are in- 
stead based on the established lists GoodTracksLoose (for charged candidates) 
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and CalorNeutral with some customized additions. The requirements for Good- 
TracksLoose, ChargedTracksAce and GoodNeutralLooseAcc are shown in Table 



Table 3.1: Criteria of some existing particle lists used in B Counting. 



GoodTracksLoosc 


ChargedTracksAce 


GoodNeutralLooseAcc 


PT > 0.05 GeV/c 
DOCA in xy < 1.5 cm 
DOCA along z < 2.5 cm 
IpI < 10 GeV/c 
Vq daughters excluded 


0.41 <0 < 2.54 


Raw energy > 0.03 GeV 
0.41 <e < 2.409 



The addition of 'Acc' (within acceptance region) to a list name denotes angu- 
lar requirements which keep the track within the detector. The list GoodTrack- 
sAccLoose is the same as GoodTracksLoosc with the additional requirement 
0.41 < 6 < 2.54 for each track. 

DOCA stands for the Distance of Closest Approach of the track to the 
interaction point. By rejecting tracks which do not closely approach the IP, a 
large number of background tracks can be removed from the event. 

Tracks clearly identified as Vq daughters are removed from the GoodTrack- 
sLoosc list. These are tracks which originate at a vertex far from the IP, and 
are obviously secondary decay products. 

3.2 Charged Candidates 

The list of charged candidates used in this thesis is based on the standard Good- 
TracksLoose with Vq daughters included. For many analyses it is convenient not 
to include them in the main charged list so they can later be identified and added 
only if desired. For example, it is beneficial for some analyses to reconstruct 
tt" mesons in this way. For B Counting, this sort of flexibility is not necessary. 
The quantities of importance for labelling events hadronic or mu-like mostly 
rely on information from the entire event such as the total event energy, the to- 
tal number of tracks, or the sphericity of the event. The sphericity is measured 
by a quantity called R2, which is the ratio of the event's second and zeroth 
Fox- Wolfram moments [TS]. 

In addition to this, some angular requirements are imposed to keep tracks 
within the detector's acceptance region. These are determined by comparing 
data and Monte-Carlo simulated data. At very high or very low angles, MC is 
less able to simulate tracks, so the maximum and minimum angular values are 
chosen to ensure data and MC agree closely in the selected region. 

At very high or very low angles there are fewer layers of the detector (espe- 
cially the drift chamber) which particles can traverse. Tracks passing near the 
edges of detector will be more difficult to reconstruct than other tracks, and 
some may be lost. This is difficult to simulate accurately with Monte Carlo. 
The angular cuts chosen are almost identical to those used in 'Acc' charged lists. 
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Figure 3.1: High and low theta regions for tracks on the GoodTracksLoose hst. 

The soUd Hne represents on-peak data, and the points (+) off-peak 
data. 



Figure 13.11 shows the high and low 9 regions for tracks on the GoodTrack- 
sLoose list. The solid curve represents on-peak data, and the dotted curve 
off-peak data. In both regions the agreement between MC and data is good 
even beyond the cuts, but sharp drop in data events at just below 0.4 and 
above 2.6 indicate the regions where the effects of reaching the detector's edges 
become apparent. The cuts are chosen on either side of these drops to ensure 
the selected region is both modelled well by MC and doesn't include regions 
where edge effects appear. 

For later convenience, this new list is given the name ChargedTracksBC. 

3.3 Neutral Candidates 

The list of neutral candidates used in this analysis is based on the standard 
CalorNeutral list. It has additional angular requirements (again to ensure the 
region studied is well-described by Monte Carlo), and a minimum energy re- 
quirement of 100 MeV. 

The reasoning behind choosing these particular angular cuts (0.42 < 9 < 
2.40) is similar to that of the charged tracks list. The high and low theta regions 
are shown in Fig. 13.21 The acceptance region for neutral particles depends on the 
angular range of the EMC (see Fig. I2.5p . so is different than for charged tracks. 
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The shape of Fig. 13.21 clearly shows the EMC construction 
are at the centre of EMC crystals. 



the highest peaks 
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Figure 3.2: High and low theta regions for candidates on the CalorNeutral list. 

The solid line represents on-peak data, and the points (+) off-peak 
data. 



The angular cut for high theta values is made to exclude the two most 
extreme crystals, ensuring neutrals are well within the detector. For the low 
theta region, the cut is based on the agreement between MC and data. Although 
EMC crystals are in place for theta values as low as 0.3, the best agreement 
between MC and data occurs above 0.42. Within reasonable limits, the exact 
choice of angular cut is arbitrary, but the range 0.42 < 6* < 2.40 is made for the 
neutral list for the above reasons. 

To reduce backgrounds, a minimum energy requirement of 100 MeV per 
cluster is also imposed on the neutral list. Erroneous neutral clusters due to 
beam gas or electronic noise have a random distribution throughout any event. 
The level of this background can be determined by examining muon-pair events 
where the muons are back-to-back. These e^e~ — > /i"*"/!" events are very clean 
— the back-to-back requirement implies the muons have not radiated photons, 
so within certain limits the muons are the only genuine physics tracks in the 
event. 

Figure 13.41 shows the distribution of energy per neutral cluster for such 
back-to-back muon events (solid curve) and compares this to mu-pair Monte 
Carlo. The events plotted pass isBCMuMu and additionally the two highest- 
momentum tracks are within 0.02 radians (1.1°) of back-to-back. A cut at 0.1 
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Figure 3.3: BB data and MC for high and low neutral angles. Candidates are 
on the Calor Neutral list. The solid line represents BB data the 
points (+) BB MC. 



GeV removes the vast majority of neutral energy from this event at the expense 
of little physics, so this requirement is included in the B Counting neutral list. 
This new list is given the name NeutralsBC. 



Tabic 3.2: Criteria of the lists used in this thesis. 



ChargedTracksBC 


NeutralsBC 


On ChargedTracks list 

PT > 0.05 GeV/c 
DOC A in xy < 1.5 cm 
DOCA along z < 2.5 cm 

IpI < 10 GcV/c 
0.42 <9 < 2.53 


On CalorNeutral list 

Cluster energy > 0.1 GeV 
0.42 <e < 2.40 
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Figure 3.4: Neutral backgrounds in Run 6 back-to-back mu-pair events. The 
solid curve represents data, and the points (+) mu-pair MC. The 
energy of each neutral cluster in an event is plotted for mu-like 
events where the muons are within 0.02 radians of back-to-back. 
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Chapter 4 



B Counting 



As its name suggests, B Counting involves determining the number of BB events 
in a particular data sample. It is necessary to know this to measure any B meson 
branching ratio for a particular final state. For example, B{B'-* D^n^) is 
approximately 0.3%. This value is equal to the number of B^ D~tt^ decays 
in a particular sample divided by the total number of B^ events in the sample. 

4.1 Method for B Counting 

The B mesons produced at BaBAR are from the decay of the T{4S) meson (a bb 
state). Due to inefficiencies, the ratio of i3°B° production to B^B^ production 
from T{4S) decay is difficult to measure, but the currently accepted value is 
consistent with 1. The BaBAR Collaboration has published the result (lOj : 

= B{T{4S) ^ B+B-)/B{T{AS) ^ BW) 

= 1.006 ± 0.036 (stat.)± 0.031 (syst.). (4.1) 

It is common to assume this number is exactly unity. Hence, to calculate a 
particular B branching ratio in a sample, it is only necessary to count the total 
number of BB events in the sample, regardless of charge. 

The current method used to count B mesons by the BaBAR Collaboration 
has been in place since 2000 [11]. In a sample of on-peak data, the number of 
BB events is equal to the total number of hadronic events {Nfj) less the number 
of non- BB hadronic events (Nx)- 

K = K^ N'^ (4.2) 

In this sense, hadronic means events which look like e~^e~ —^ hadrons. The 
superscript "0" in each case indicates these quantities are the actual numbers 
produced of each event type. In reality, every quantity has an associated effi- 
ciency since nothing can be counted perfectly. In the B Counting code, we can 
only use the number of each quantity counted, and the number of B mesons 
produced is calculated from the number of B mesons counted by dividing by 
the BB efficiency Eb- 

The number of on-peak continuum events can be found by scaling (by lu- 
minosity) an off-peak sample where, without T(45') production, all events are 
non-BB. The small decrease in energy from on-peak to off-peak data-taking 
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changes all continuum production rates slightly, but almost all of these events 
scale in the same way as muon-pair events (i.e. e+e" jX^ Hence, the 
ratio of muon-pairs in the on- and off-peak samples is approximately equal to 
the ratio of the luminosities of the two samples. B Counting relies on the fact 
that the ratio of muon-pair to continuum (i.e. non-BB) events in a sample is 
the same (to an excellent approximation) regardless of the CM energy: 

In this thesis, unless stated otherwise, primed symbols represent off-peak values. 



4.2 Derivation of the B Counting Formula 

Suppose we wish to calculate the number of B mesons in a sample of on-peak 
data of luminosity C, using an off-peak data sample of luminosity C . For sim- 
plicity, off-peak quantities are primed, and hadronic, mu-pair, BB and contin- 
uum quantities have the subscripts H, /x, B and X respectively. The numbers 
produced of each quantity (as opposed to those counted by the B Counting 
selectors) have the superscript "0" . 

In a sample of on-peak data, the number of BB events is equal to the total 
number of hadronic events {Njj) less the number of non-BB hadronic events 

K=N°j- N°x (4.4) 

The number of on-peak continuum events can be found by scaling (by lu- 
minosity) an off-peak sample (in which all events arc non-_Bi3). The small 
decrease in energy from on-peak to off-peak data-taking changes all continuum 
production rates slightly, but almost all of these events scale in similar ways 
with luminosity. 

For any particular type of event, the number counted is equal to the number 
produced multiplied by the efficiency (e). So for example, 

and 

%=e'^%'=sy^C'. (4.6) 

For off-peak data, there is no Y{4S) production, so we can assume that all 
hadrons are from continuum events (and hence the symbols N'jj and N'x are 
equivalent): 

K=N'x= s'xN'° = e'xa'xC. (4.7) 



We define: 



SuCTfj. ^X^X 
= K^-Kx. (4.9) 
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We combine (g^, (jiT)) and (g^ to give: 

#.= ^^.= ^. (4.10) 

The hadronic events in the on-peak sample consist of continuum and BB 
events: 

iV^=7V^+iV^ (4.11) 
and the number of counted hadronic events is 

Nh = exN"x+SBN"B- (4-12) 
Hence from (|4.10p . we can write 

ebNI, = Nh-exN^'x (4.13) 
= NH-ex(yxC (4.14) 



N,, 



= Nn-'-^-^-e^a^--^ (4.15) 



N'rr 

= Nh-N^--^-k. (4.16) 

Hence, the number of BB mesons produced in the on-peak sample is given 

by: 

N% = —{NH-N^-Roff-K), (4.17) 

where 

AT'' 

Rob ^ (4.18) 

In general, k is close to unity. The exact values and uncertainties of and k,x 
are discussed in detail later. 



4.3 B Counting Goals 

Figure 14.11 demonstrates the scaling of on- and off-peak data by the number of 
muon-pairs in each sample. The figure is a plot of total energy for events passing 
all hadronic cuts, except for one on ETotal itself (the existing B Counting cut is 
set at 4.5 GeV). Samples of on- and off-peak data from the same time period are 
plotted — in this case from 'Run 4', which took place between September 2003 
and July 2004. The off-peak data is luminosity-scaled (by the ratio of muon 
pairs in each sample) to match the on-peak data. 

When correctly-scaled off-peak data is subtracted from on-peak data, the 
remaining events are those from T{AS) decay, the vast majority of which are 
BB events. This BB data is shown in Fig. 14. II and again in Fig. 14.21 where it is 
overlaid with BB Monte Carlo (MC) generated events. The agreement between 
data and MC is generally very good. 

Ideally the B Counting hadronic selection efficiency should be: 
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Figure 4.1: Total energy divided by CM energy (ETotal/eCM) 
for events passing all B Counting hadronic cuts apart from one on ETotal for 
a sample of data from Runs 5 and 6. The bin width is 0.007. The solid (+) 
histogram represents on-peak (off-peak) data. The off-peak data is scaled by lu- 
minosity to the on-peak data, and the dotted histogram is the difference between 
the two. 
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Figure 4.2: ETotal/cCM on- and off-peak subtraction, overlaid with BB Monte 
Carlo. The solid histogram represents data (off-peak subtracted 
from on-peak in Fig. 14.11) and the dotted histogram is from MC. 
The histograms are normalised to have the same area from 0.0-1.0. 
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• high for B meson events, 

• insensitive to changing background conditions (i.e. stable in time), and 

• in good agreement with simulation BB Monte Carlo. 

Similarly the B Counting muon selection should also be insensitive to back- 
grounds and be in good agreement with MC. The specific selection requirements 
are described in later chapters. 

More than 400 million BB pairs have been created in BaBAR during its 
running, so B Counting is a very high statistics task and uncertainties are 
predominantly systematic. The main aim of the work described in this thesis is 
to find selection criteria optimising these requirements and minimising the final 
systematic error in the number of B pairs counted. 

4.4 TrkFixup 

This thesis proposes the first major update to B Counting since the original 
method was devised in 2000. Since then, many improvements have been made 

to both hardware and software within BaBAR which have improved the precision 
of many analyses. One of the most significant for B Counting, named TrkFixup 
was proposed in 2005 and fully implemented by mid-2007. 

A systematic uncertainty present in any BaBAR analysis comes from the 
tracking precision. There are many reasons the tracking information of charged 
particles passing through the detector can be incorrect. For example, a collection 
of drift chamber and SVT hits due to noise could be incorrectly reconstructed as 
a charged track. Similarly low-momentum particles which loop around within 
the detector can be misidentified as more than one track. 

TrkFixup improves the tracking precision through several methods. Firstly, 
it improves track resolution by modifying pattern recognition algorithms. It also 
increases the track finding probability and removes unwanted ('junk') tracks 
such as duplicate (e.g. looping) tracks and backgrounds. This is achieved by 
re-examining 'questionable' tracks with more sophisticated algorithms than pre- 
viously used. 

These improvements arc important to B Counting for several reasons. Re- 
ducing sensitivity to beam backgrounds is a key part of this thesis and the 
addition of TrkFixup to track reconstruction improves this even before any ad- 
justment to the B Counting cuts is made. Some other improvements are obvious, 
such as increased precision in measurements of track momentum and position. 
This is important for border-line events which narrowly pass or fail a particular 
cut. 

When counting muon-like events, a requirement of the B Counting code is 
that the event have at least two 'good' tracks. TrkFixup affects B Counting 
directly through this cut, as events previously with only one good track could 
now be found to have two or more. Similarly, an event previously labelled as 
mu-like could be found after TrkFixup to have fewer than two good tracks and 
thus fail the selection. 
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4.5 Original B Counting Requirements 

The original requirements (i.e. those in place before the research described 
herein) for an event to be labelled as mu-like or hadronic for B Counting 
purposes were proposed in 2000 [11] . An event passing these criteria during 
processing receives a label (known as a 'tag') isBCMuMu or isBCMultiHadron 
respectively. The criteria are summarized in Table 14.11 

During reconstruction, events can receive one or more BGF tag. They have 
very simple criteria and are designed to quickly accept all possibly useful physics 
events. The BGF tags enable some time to be saved during reconstruction, as 
an event must have at least one of them to be fully processed. BGFMultiHadron 
and BGFMuMu are two of these which are used in B Counting. 

BGFMultiHadron has the very simple selection requirements of at least three 
charged tracks, and R2 of less than 0.98 in the CM frame. R2 is the second 
Fox- Wolfram Moment divided by the zeroth Fox- Wolfram Moment and is a 
measure of the event's sphericity [T3]. A pair of back-to-back particles have 
R2 approaching one, while very spherical events (i.e. ones with many tracks 
distributed evenly around the event) have R2 approaching zero. Thus the events 
which do not receive the BGFMultiHadron tag are mostly e+e~ — > t^t^ . The 
three track requirement also filters out several types of two-photon and Initial 
State Radiation (ISR) events. 



Table 4.1: Criteria of the isBCMuMu and isBCMultiHadron B Counting tags. 



The definitions of each quantity are given in the text. 



isBCMuMu Tag 


isBCMultiHadron Tag 


BGFMuMu 
masspair > 7.5 GeV/c^ 
acolincm < 0.17453 (10°) 

nTracks > 2 
maxpCosTheta < 0.7485 
maxpEmcCandEnergy > GeV 
maxpEmcLab < 1 GeV 


BGFMultiHadron 
R2AU < 0.5 
nCTL > 3 
ETotal > 4.5 GeV 
PrimVtxdr < 0.5 cm 
|PrimVtxdz| < 6.0 cm 



The criteria of BGFMuMu are also very broad, and are designed to quickly 
reject events which cannot be mu-pair events. To be labelled BGFMuMu, the 
two tracks of highest momenta must have momenta above 4 GeV/c and 2 GeV/c 
respectively and be approximately back-to-back (2.8 < 0i + 6*2 < 3.5). Also, the 
energy deposited in the EMC by the two highest-momentum tracks must add 
to less than 2 GeV. 

The isBCMuMu tag has many criteria which depend on the properties of the 
two highest-momenta tracks — if the event is a genuine mu-pair event, these will 
almost certainly be the two muons. The quantity masspair is the invariant mass 
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of the pair and acolincm is the acohnearity (the angular departure from back- 
to-back) of the pair in the centre-of-mass (CM) frame. Also in the CM frame, 
both tracks must have | cos(0)| < 0.7485. The maximum energy deposited in the 
EMC by one of these tracks must be less than 1 GeV in the laboratory frame, 
and at least one must leave some energy in the EMC. The former condition 
rejects approximately 99.7% of all Bhabha events, which on average deposit 
just over 5 GeV of energy per electron track. The final isBCMuMu quantity is 
nTracks, the number of reconstructed charged tracks on the TaggingList list. 

The criteria of isBCMultiHadron mostly involve quantities related to the 
entire event. R2AII is R2 calculated only using charged and neutral candidates 
from within the detector acceptance region (i.e. from the ChargedTracksAcc 
and GoodNeutralLooseAcc lists). There are also requirements on the number of 
tracks on the GoodTracksLoose list (nGTL) and the total event energy recorded 
in the laboratory frame (ETotal). 

A possible source of background events is beam-wall interactions. These are 
showers caused by a particle from the beam interacting with the walls of the 
beam pipe and can resemble genuine physics events. The primary vertex is the 
point calculated after reconstruction to be the most likely vertex for all tracks in 
the event. To reject beam-gas interactions, isBCMultiHadron has requirements 
on the distance between the primary vertex and the beamspot in the xy plane 
(PrimVtxdr) and in z (PrimVtxdz) . 

The number of B mesons in a particular on-peak sample can be calculated 
using these tags and (|4.17p by equating the number of isBCMuMu events with 
Nf^ and the number of isBCMultiHadron events with Nh- 

The value of i?off depends on beam conditions and backgrounds and so 
changes with time. Off-peak running occurs only for short periods around larger 
blocks of on-peak running, so when calculating the number of B mesons in any 
particular on-peak sample, the value of Rojf is based on the closest blocks of 
off-peak running. 

Many effects including changing background conditions cause long-term time- 
variations in selector (hadronic and muon-pair) efficiencies and this is discussed 
in the following chapters. The association of on-peak data with only the most 
recent off-peak running period ensures this variation has a minimal effect on B 
Counting. 
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The motivation for the work in this thesis is to improve the understanding of 
B Counting at BaBAR and if possible to modify the method to reduce overall 
uncertainty. With more than 400 million BB pairs recorded, statistical effects 
are often negligible, so improvements must come through reducing systematic 
uncertainty. This chapter describes proposed modifications to the mu-pair part 
of the B Counting code. 

5.1 Monte Carlo Simulation 

Muon-pair events are simulated with the KK2F MC generator . Before any 
detector simulation is applied, the physics cross-sections for the process e^e~ 
fX^ within the detector can be calculated using a BABAR application called 
GeneratorsQAApp. Based on 10 million generated on- and off-peak events, the 
cross-sections are (1.11853±0.00011) nb and (1.12647±0.00011) nb respectively. 

For detector simulation of all MC events, the BaBAR Collaboration uses 
GEANT 4 [H]. After GEANT has been applied to the simulated events, they are 
made to correspond to a particular time by adding background frames from data. 
These frames are random 'snapshots' of detector occupancy and are recorded 
at a rate of 1 Hz. The events are then reconstructed in the same way as data. 
This background mixing increases the realism of the Monte Carlo simulation 
by more accurately reflecting the processes occurring simultaneously within the 
detector. 

The detector conditions (such as the number of non-functioning EMC crys- 
tals), the levels of beam-background, luminosity and other effects vary with 
time and cause a time-dependence in the efficiency of the mu-pair selector. The 
efficiency of isBCMuMu (and almost any other quantity) for a sample of MC 
depends on the time at which it was 'taken', i.e. on the time of the mixed-in 
backgrounds. Here, efficiency means the proportion of generated mu-pair events 
which pass the mu-pair B Counting selector. In this way MC simulates the same 
effect in data. 

The efficiency of isBCMuMu for simulated mu-pair MC is approximately 
43%. As described in Chapter |4l the number of mu-pair events in a data sample 
is used to scale the sample by luminosity. In modifying the muon selection, the 
goal is to reduce the total uncertainty when counting B mesons, so the absolute 
value of the efficiency is not as important its variation in time. 
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5.2 Cut on Muon Pair Invariant Mass 

Ideally the muon selection will function identically on both on- and off-peak 
data. The total change in energy between the two is small (around 40 MeV), 
but any cut which varies depending on CM energy should be altered to minimize 
the effect. 

Of the isBCMuMu cuts listed in Table WA\ the one altered the most by this 
effect is masspair — the invariant mass of the two highest momentum tracks. 
This cut is designed to reject tau pairs but accept almost all genuine muon 
events. The tau pairs are very short-lived and decay into multiple particles 
before reaching the beam-pipe walls. These events fail because the two highest- 
momentum tracks do not carry the majority of the energy and their invariant 
mass is comparatively low. For muon-pair events, masspair peaks near the beam 
energy. 

A simple method of removing sensitivity to CM energy (eCM) is to define 
a new quantity, masspair scaled by the total energy — masspair/eCM — and 
cut instead on this. If the cut value in Table BT] is scaled by the on-peak beam 
energy of 10.58 GeV, the revised cut is: 

masspair 7.5 GeV/c^ 

-^CM- > 10.58 GeV " ("'^ 
In this thesis, we set c — 1 for convenience unless otherwise stated. 



5.3 Minimum EMC Energy Requirement 

As described in Chapter 01 to receive the isBCMuMu tag an event must meet 
certain EMC requirements. Of the two highest-momenta tracks, the maximum 
energy deposited in the EMC by either one must be non-zero and less than 1 
GeV. It was found during this research that depending on the time-of-running, 
up to 2% of events failed muon selection because both tracks deposited no energy 
in the EMC, i.e. events which would otherwise pass the selection were rejected 
because the quantity maxpEmcCandEnergy was zero. 

In the original B Counting code, the maximum EMC energy cut was calcu- 
lated in the lab frame, while the minimum was calculated in the CM. However, 
any zero measurement in the EMC is equivalent regardless of frame, so for sim- 
plicity in this thesis, the minimum and maximum allowed EMC deposits are 
both calculated in the lab frame. 

This effect was studied in samples of at least a million on-peak data events 
with the BGFMuMu tag from each of the six major running periods (Runs 1- 
6). Some of the quantities of interest for events passing all other muon selection 
requirements are shown in Tables 15.11 and 15.21 The value of plEmcCandEnergy 
is the amount of energy left in the EMC by the highest momentum track, and 
p2EmcCandEnergy the amount for the second highest momentum track. The 
maximum of the two is named maxpEmcCandEnergy. The distribution of this 
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Figure 5.1: Distribution of uiaxpEmcCandEnergy for on-peak Run 6 data pass- 
ing all other muon-pair selection cuts. 



quantity for a sample of Run 6 on-peak data passing all other muon-selection 
cuts is shown in Fig. 15.11 

Any track has a non-zero chance of leaving no energy in the EMC. This could 
be due to any of several reasons including reconstruction errors or the particle 
hitting a gap between two EMC crystals. Assuming the individual probabilities 
are independent, the probability of both of the two highest momentum tracks 
leaving no EMC energy should be approximately the square of this number. 
However, as can be seen in Tables \5A] and [5T2| the number of events with both 
EMC clusters equal to zero (i.e. where maxpEmcCandEnergy=0) is much higher 
than would be expected. 

This effect does not appear to be due to EMC failures in a particular region. 
Figure 15.21 shows the angular distribution in cos{9) and of the muons in a 
sample of Run 6 data events which pass all other muon-selection cuts, but 
which have maxpEmcCandEnergy of zero. Note that the distribution is fairly 
randomly spread in 9 and (f>, although some broad bands are visible at constant 
(j) due to detector geometry. 

The high proportion of events with maxpEmcCandEnergy of zero indicates 
some other effect is present — it is especially noticeable in Runs 2 and 3 where 
the probability of both tracks leaving nothing in the EMC is greater than the 
probability that the highest momentum track alone deposits nothing. This was 
examined by Rainer Bartoldus and the BaBAR Trigger, Filter and Luminosity 
Analysis Working Group (AWG). The predominant cause was found in events 
where the DCH and EMC event times differed. 

The Drift Chamber and EMC have different timing systems, so can occa- 
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Figure 5.2: Angular distribution of mucus from Run 6 data with maxpEmc- 
CandEnergy of zero. The value of cj) and cos(6') in the CM frame 
are plotted for the two highest-momentum tracks in events which 
pass all other muon-selection requirements. The number of events 
in each bin of this 100 x 100 grid is proportional to the area of the 
box at the bin-centre. 



sionally disagree on the event time (tg). The EMC tQ is calculated in two stages. 
First, the average arrival-time of energy in a wide time- window around the event 
is found. The window is then narrowed around this average to select the EMC 
hits to be reconstructed. This method is susceptible to calculating incorrect 
^0 if background events such as stale Bhabhas arrive at an inopportune time. 
Stale Bhabhas are unwanted Bhabha events which happen to occur at almost 
the same time as another physics event. 

Such an event passing through the EMC at either end of the event window 
could cause the average arrival time of energy to shift to a time when no energy 
is in the EMC. Through this process events that would otherwise pass muon 
selection can have no energy visible in the EMC and be rejected. The proportion 
of events with EMC information lost in this way is complicated to simulate, but 
Tables 15.11 and 15.21 show reasonable agreement (usually within a factor of 2) 
between data and Monte Carlo. 

The majority of these events can be regained by changing the minimum 
maxpEmcCandEnergy to: 

• Either maxpEmcCandEnergy > 0, 

• or at least one of the two highest momentum tracks is identified as a muon 
in the IFR. 



Table 5.1: Effect of the minimum maxpEmcCandEnergy requirement on samples from Runs 1-3. 





Run 1 




Run 2 




Run 3 






Data 


MC 


Data 


MC 


Data 


MC 


# BGFMuMu Events 


1.001 X lO'' 


0.807 X lO*" 


1.542 X lO** 


2.580 X lO** 


4.622 X lO'' 


1.718 X 10^ 


rass all cuts but minimum 


U.Dol X 10 


0.576 X 10 


1.055 X 10 


1.837 X 10 


3.155 X 10 


1.223 X 10 


EMC energy (f) 














% of t with 


0.52 


0.32 


1.34 


0.41 


1.09 


0.38 


plEmcCandEnergy=0 














and p2EmcCandEnergy^ 














% of t with 


0.91 


0.70 


1.71 


0.83 


1.44 


0.79 


p2EmcCandEnergy=0 














and plEmcCandEnergy^ 














% of t with 


0.04 


0.14 


1.70 


1.91 


1.37 


1.66 


maxpEmcCandEnergy=0 (|) 














% of 1 with at 


23.3 


85.4 


83.3 


93.4 


80.6 


91.3 


least one track identified 














as a muon in IFR 















Table 5.2: Effect of the minimum maxpEmcCandEnergy requirement on samples from Runs 4-6. 





Run 4 




Run 5 




Run 6 






Data 


MC 


Data 


MC 


Data 


MC 


# BGFMuMu Events 


2.516 X 10'' 


3.556 X 10^ 


2.694 X 10'' 


3.045 X lO** 


3.151 X 10*' 


1.003 X lO** 


rass all cuts but minimum 


1.724 X ID 


2.542 X 10 


1.874 X 10 


2.186 X 10 


2.154 X 10 


A Tl 1 N , 1 a6 

0.711 X 10 


EMC energy (f) 














% of t with 


0.66 


0.45 


0.65 


0.43 


0.63 


0.36 


plEmcCandEnergy=0 














and p2EmcCandEnergy^ 














% of t with 


1.07 


0.84 


1.08 


0.83 


1.04 


0.76 


p2EmcCandEnergy=0 














and plEmcCandEnergy^ 














% of t with 


0.19 


0.27 


0.16 


0.31 


0.21 


0.42 


maxpEmcCandEnergy=0 (|) 














% of 1 with at 


65.1 


93.4 


89.8 


97.9 


87.6 


93.5 


least one track identified 














as a muon in IFR 
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A requirement on the IFR has not been used previously in B Counting. This 
is because in general its muon identification efficiency is too low and variability 
with time too large to provide stable identification. However, in this case the 
IFR is used only to recover genuine muon-pair events that would otherwise be 
lost, so enables a net efficiency gain and improved stability in time. 

Muons are identified in the IFR by a Neutral Net muon selector [TH] , specifi- 
cally by muNNVeryLoose. The efficiency of this selector varies with time, but in 
muon-pair MC the probability that at least one of the two highest-momentum 
tracks passes is typically over 90%. 

The final row of Tables [01 and [5T^ shows the percentage of events recovered 
by adding this IFR requirement. In Run 6 data for example, 87.6% of events 
which pass all other isBCMuMu requirements but which have maxpEmcCan- 
dEnergy of zero are recovered. 

5.4 Optimization of Cuts 

The muons counted by isBCMuMu selection are used to measure the luminosity 
of a sample of events. Since on- and off-peak data are taken at different times, 
the two have different background conditions. Hence, to correctly perform the 
B Counting background subtraction it is important the muon-pair selection be 
as insensitive as possible to changing backgrounds. The variation is reduced 
to some extent by changing to the particle lists described in Chapter [3] and by 
altering the masspair and maxpEmcCandEnergy cuts as described above. 

This sensitivity can be measured by the changes in muon-pair MC efficiency 
with time. Ideally the efficiency should be identical regardless of run number — 
a sample of generated mu-pair events with Run 1 backgrounds should have the 
same efficiency as the same events with Run 6 backgrounds. In practice this is 
not the case, but by choosing different cut values for the requirements outlined 
in Table |4T] it can be reduced. 

The requirement of BGFMuMu is essential as it quickly filters out large 
amounts of events which cannot be muon-pairs. Likewise there is little flexibility 
in the requirement of at least two well-identified charged tracks. The selector 
should identify events with well-tracked muons, and all events of this type will 
have at least two tracks on the charged list. The modification of the minimum 
maxpEmcCandEnergy requirement described above helps veto Bhabha events 
and is not further adjusted. 

The four remaining isBCMuMu requirements have slightly more fiexibility 
in their cut values. The justification for the current cuts is given in [IT, but the 
addition of TrkFixup and the different particle lists (ChargedTracksBC and Neu- 
tralsBC) used in this thesis are reason enough that they should be re-evaluated 
at this stage. Additionally, the method used in this thesis to find an optimal 
cut-set is improved and more robust. 

To optimise the isBCMuMu cutset, the month-by-month efficiency for muon- 
pair MC in both on- and off-peak running was calculated for different cutset 
choices. Up to and including Run 6, there are 63 months with on-peak muon- 
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pair MC and 27 months with off-peak MC. For convenience, each month from 
the time of running is given an integer label — from January 1999 (month 0) 
to December 2007 (month 107). 

The total number of cutsets trialled in the optimization procedure was 13 x 
17x 20 X 14 = 61880. The weighted mean efficiencies of these cutsets are between 
0.262 and 0.450. 



Table 5.3: Cut optimization of isBCMuMu. 



Quantity 


No. of increments 


Range 


masspair/eCM 


13 


(6.0 - 9.0)/10.58 


maxpCosTheta 


17 


0.650 - 0.765 


acolincm 


20 


0.01 - 0.20 


maxpEmcLab 


14 


0.7-2.0 GeV 



The forward and backward ends of the detector have different acceptance 
regions, and the different angular cuts on ChargedTracksBC and NeutralsBC 
reflect this. When varying maxpCosTheta (the maximum value of | cos 9cai\ for 
the two highest momentum tracks), some care must be taken to ensure the cut 
is not varied into a region excluded by the list's own angular cuts. 

To find the 9cm angles corresponding to the angular cuts of ChargedTracksBC 
for a typical muon we boost into the laboratory frame from the CM by f3j = 
—0.56 along z. The method for calculating the transformation of the angles of 
a velocity can be found in many relativistic kinematics texts, for example |20| . 
We find: 

vsinO ^. 

tanOcM = —, -3 , ^ ^ 5.2) 

7(1; cos^ + pc) 

where v is the velocity of a particle travelling at angle 6 relative to the z axis 
in the laboratory frame. Practically all muons from e+e^ /^^M^ have CM 
momenta between 5.0 and 5.5 GeV/c, which corresponds to having velocities 
of 0.9999c in the laboratory frame (the two are identical to four significant 
figures). Using (|5.2p and this velocity, the minimum and maximum angles of 
ChargedTracksBC (0.42 and 2.53) boost to 0.698 (cos ^cm = 0.766) and 2.776 
(cos 6'cM = 0.934) in the laboratory frame respectively. 

Hence, when varying maxpCosTheta, values above 0.766 should be ignored 
to agree with the angular requirements of the charged list. In practice, cut 
values up to and including 0.765 are trialled. 

The on- and off-peak muon-pair MC efficiency variation is shown in Fig. 15.31 
and l5.4l The cutset shown has fairly typical time- variation. It has an acolincm 
cut of 0.07, masspair/eCM of (8.00/10.58), maxpCosTheta of 0.765 and a 
maximum maxpEmcCandEnergy cut of 1.0 GeV. 

In each month i, the statistical error A(£i) is assumed to be binomial with 



(5.3) 
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Figure 5.3: On-peak mu-pair MC efficiency variation. Details of cutset choice 
and error bars are in text. The value is 124.04 and M (the 
number of degrees of freedom) is 63. The dashed lines show ±0.5% 
from the mean. 
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Figure 5.4: Off-peak mu-pair MC efficiency variation for the same cutset choice 
as on-peak distribution. The value is 47.94 and M (the number 
of degrees of freedom) is 27. The dashed lines show ±0.5% from the 
mean. 
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where Si is the efficiency and Ni the total number of generated events for that 
month. If we assume the efficiencies come from a statistical distribution which 
is constant in time, the weighted mean e (the solid horizontal line on the plots) 
is given by 

1 

i=l 

where M is the number of months with MC events, Wi = l/(A(£i))^ and w = 
J2 Wi ■ The statistical standard deviation (shown by dashed horizontal lines) 
of e is The variation around the weighted mean is obviously not 

random, and this shows the presence of some non-statistical systematic effects. 

The amount of systematic variation present can be measured by calculating 
the statistic 

for each cutset. If the variation was entirely statistical, an average value of 
M — 1 would be expected. The value of M — 1 is 62 for on-peak MC and 26 for 
off-peak MC. For the cutsets trialled, is much higher than this: for example 
124.04 and 47.94 in Fig. [Ol and EUl The on- and off-peak distributions of x^ 
values for the 61880 cutsets are shown in Fig. 15.51 and l5. 61 

All cutsets trialled have a lower x^ statistic than the existing B Counting 
muon-pair selection. For the same samples of on- and off-peak MC, the x^ 
statistics of the isBCMuMu tag were 1016.5 and 402.25 respectively. 



5.5 Estimating Systematic Uncertainty 

From Fig. 15.51 and 15.61 it is clear that all cutsets have values indicating vari- 
ation beyond statistical fluctuations. There is however, a broad range of x^ 
values — this indicates that some choices of cutset are objectively better (i.e. 
less sensitive to systematic effects) than others. 

To gauge the size of the systematic uncertainty we look for cIq in a modified 

statistic 

such that 

x''(do) = M-l. (5.7) 

Equation (j5.6p is equivalent to a polynomial of order 2M in d. Clearly for 
d = 0, x'^(O) ^ > M -1 and for large d, x'^{d) ^ 0. So, since jSS]) is 
continuous, by the Intermediate Value Theorem there is at least one real do 
satisfying (|5.7p . 

In practice, do is found numerically by incrementing d from zero in steps of 
1 X 10^^ until x''^{d) < M — 1. The value of d at this point is taken to be equal 
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Chi-squared statistic 



Figure 5.5: Distribution of statistic for trialled cutsets with on-peak mu-pair 
MC. M (the number of degrees of freedom) is 63. 




Chi-squared statistic 



Figure 5.6: Distribution of statistic for trialled cutsets with off-peak mu-pair 
MC. M (the number of degrees of freedom) is 27. 
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to do- The ranges of do values for different cutsets are (5.57 - 9.27) x 10 in 
on-peak MC and (2.93 ~ 7.89) x 10"'' in off-peak MC. 

The total uncertainty for a particular cut-set is then the sum in quadrature 
of statistical uncertainty (l/y^) and systematic uncertainty (do). The 'optimal' 
set of cuts is the one which causes the lowest uncertainty in the total number 
of B mesons. 

It is important to note the difference between the uncertainty used for the 
purposes of optimization and the value of the systematic uncertainty in data. 
The optimization procedure seeks to find the cutsets (muon-pair and hadronic) 
which cause the lowest uncertainty in B Counting due to time- variation of MC. 
This amounts to finding the cutsets which are most stable in time. 

The systematic uncertainty for MC does not correspond directly to the same 
uncertainty in data. Steps are taken in data to reduce the effect of varying effi- 
ciencies, and the association of on-peak data with only the most recent Roff value 
reduces this effect to a large degree. The estimation of systematic uncertainties 
in data is discussed in Chapter [9l 



5.6 Variation in 

The results above can be used to evaluate the value of as defined in (|4.8p . In 
past B Counting code, k, the product of and kx has been set equal to unity. 
Kfj, is equal to the ratio of off- to on-peak muon-pair cross-sections, multiplied 
by the respective efficiencies. 

The theoretical cross-section for electron annihilation to a muon-pair e^e~ —^ 
at CM energy squared of s is approximately [22]: 

47ra2 86.8 nb 
OS s (m Gev ) 

So, since the muon-pair cross-section scales inversely with s, the expected value 
of is approximately: 

= ^ - wS - 1-0076. (5.9) 
e^cr^ 10.54^ 

To find the time-variation in some care needs to be taken when choosing 
time bins. On- and off-peak running occur at different times, and often the gap 
between off-peak running periods is quite large. To manage this, the running 
time is divided into 24 unequal bins such that each contains periods of both off- 
and on-peak running. This makes it possible to calculate a value for in each 
time bin. 

Variation in is calculated using the same method used to calculate varia- 
tions in mu-pair efficiency. The statistical uncertainty in includes the effects 
of the independent statistical uncertainties of on- and off-peak mu-pair MC ef- 
ficiencies which are given by (j5.3[) . There is also a small contribution from the 
uncertainty in the MC cross-sections, A((t^) and A((T^). 
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Figure 5.7: Time-variation of for the cutset shown in Fig. 15.31 and 15.41 The 
value is 21.01, M is 24 and the weighted average is 1.00854. The 
dashed hues show ±0.5% from the mean. 



(5.10) 



and similarly for off-peak MC. The method used for calculating propagation 
of uncertainties is given in Appendix [Xj The statistical uncertainty in for 
time-period i is then given by (|A.4p : 



A(K^,i) 



ey^ ■ A(£,cr^) 



(5.11) 



The shape of the time distribution for the same cutset considered in Section 
is shown in Fig. 15.71 It is much closer to a constant than for either the on- 
or off-peak efficiencies, and the variation in time seems to be mostly consistent 
with only statistical variation. This is an indication that much of the variation 
in mu-pair MC is due to systematic effects which affect on-peak and off-peak 
MC similarly. 

With 24 periods of time considered, the value of Af — 1 for is 23. The 
statistic can again be calculated, and for the pictured cutset it has a value of 
21.01. The weighted mean value of k,^ for this cutset is 



= 1.00854 ±0.00044 



(5.12) 
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which within two (statistical) standard deviations of the theoretical expectation. 
The distribution of values for all cutsets is shown in Fig. 15.81 Across all 
cutsets, the mean value of is 24.88. 




Chi-squared statistic 



Figure 5.8: Distribution of statistic for all trialled cutsets. M (the number 
of degrees of freedom) is 24. 

As before, the amount of systematic uncertainty present in each cutset's 
variation can be estimated by finding do satisfying a modified (|5.7p . We do not 
allow (P to be negative, so for cutsets with less than 23, c?o is set to zero. 
Following the same method as above, the range of do values is between zero and 
0.00172. This is added in quadrature with the statistical uncertainty to give the 
total uncertainty on each cutset's value. 
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The criteria for good hadronic selection requirements were mentioned in Section 
14.11 This chapter describes proposed modifications to the hadronic part of the 
B Counting code. 

The change to the NeutralsBC fist with a minimum neutral energy require- 
ment reduces the sensitivity to backgrounds before any optimization is per- 
formed on the hadronic selection. Several of the other cut quantities are changed 
to match the new lists. Instead of requiring at least three tracks on the Good- 
TracksLoose list, the cut is instead based on ChargedTracksBC. Previously, the 
values of ETotal and R2A11 were based on the the GoodNeutralLooseAcc and 
ChargedTracksAcc lists. Now, they are calculated by combining all candidates 
on the ChargedTracksBC and NeutralsBC lists. 

Similarly, PrimVtxdr and PrimVtxdz (the distance of the event's primary 
vertex from the beamspot in r and z) were formerly calculated by finding the 
vertex of all tracks on ChargedTracksAcc. In this thesis, these quantities are 
calculated from tracks on an intermediate charged list, which is equivalent to 
ChargedTracksBC before any cuts on DOCA (in either xy or z) are made. 

6.1 Monte Carlo Simulation 

The method for optimizing hadronic selection is similar to the mu-pair selection 
optimization discussed in the previous chapter. Hadronic continuum events 
are simulated with e'^e~ — s- uu/dd/ss and e^e~ — > cc on-peak and off-peak 
MC using the JETSET MC generator [23]. As with muon-pair MC, once the 
MC events are generated, detector response is simulated with GEANT 4 and 
background frames are added. 

In data, there are additional sources of hadronic continuum events — Initial 
State Radiation (ISR) and two-photon events. These are small in comparison 
to the other continuum events and their effects are discussed in later chapters. 
The present optimization procedure seeks to minimize the time-sensitivity of 
hadronic selection for uds and cc Monte Carlo only. 

The production cross-sections for the continuum MC used in the optimiza- 
tion procedure are given in Table 16.11 These values are based on 1 million 
on- and off-peak events of each type simulated by the KK2F generator. BB 
events are simulated in a similar way to continuum MC, but with the EvtGen 
MC generator [21]. The efficiency of the isBCMuhiHadron tag for BB MC is 
approximately 96%. 
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Table 6.1: Continuum MC cross-sections used in optimization. 



Continuum 
MC Type 


On-Peak a 


Off-Pcak cr 


uu 


1.5735 ±0.0013 nb 


1.5847 ±0.0013 nb 


dd 


0.3931 ± 0.0003 nb 


0.3959 ±0.0003 nb 


ss 


0.3751 ± 0.0003 nb 


0.3778 ±0.0003 nb 


cc 


1.2929 ± 0.0004 nb 


1.3021 ± 0.0004 nb 



6.2 Cut on Total Energy 

As in Section [521 hadronic selection should function in the same way, regardless 
of CM energy. Of the isBCMultiHadron cuts listed in Table the one affected 
by this is the cut on total event energy in the laboratory frame. This dependence 
is removed by scaling by eCM. If the cut value is scaled by the on-peak beam 
energy, the adjusted requirement is: 

ETotal 4.5 GeV „ ^ ,^ , 

The cut on total lab energy is primarily in place to reject beam-gas and 
two-photon events. The adjusted cut slightly increases the isBCMultiHadron 
efRciency for off-peak continuum events. 



6.3 Cut on Highest-Momentum Track 

Events of the type e~'"e~ £~^£~ where £ is either e or n typically have at least 
one track with approximately half the beam energy in the CM frame. This does 
not happen for B decays, as the B mesons very quickly decay into other particles 
well within the detector. Typically each B meson decays into two tracks, each 
with approximately a quarter of the beam energy. One possible way to increase 
the Bhabha and mu-pair rejection in the hadronic selection without affecting 
BB efficiency is to reject events where the highest momentum track has close 
to half the beam energy. 

Some poorly-tracked particles have an artificially high momentum (for ex- 
ample a momentum of 4 ± 2 GeV/c), and SVT-only tracks are especially prone 
to this. At short distances chord and arc length are very similar, so to precisely 
measure a track's curvature (and thus momentum) it is usually necessary to 
make measurements outside the SVT, some distance from the beamspot. A 
minimum requirement for this is at least one hit in the drift chamber. 

In hadronic selection optimization, one variable is a new proposed cut on 
plMag, the CM momentum of the highest-momentum track. If the highest- 
momentum track has at least one hit in the drift chamber and its momentum 
is greater than the cut value, the event is rejected. 

Figure 16.11 shows the distribution of plMag for BB MC events where the 
highest-momentum track has at least one drift chamber hit. The sample contains 
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1.8 million Run 6 MC events. Clearly, above 3.0 GeV/c the number of events 
is negligible. 




Figure 6.1: Distribution of plMag for a sample of BB Monte Carlo events where 
the highest-momentum track has at least one DCH hit. 



In the optimization procedure, the cut on plMag (with the requirement on 
DCH hits) is varied between 3.0 and 5.0 GeV/c or is turned off altogether. 

6.4 Optimization of Cuts 

Similarly to mu-pair selection, there are certain cuts necessary for isBCMul- 
tiHadron which allow little flexibility. These include the cuts on PrimVtxdr 
and PrimVtxdz. The reasons for the current cut values on these quantities are 
given in [iT]. Both cuts are designed to reject backgrounds, while maintaining 
high BB efficiency. They especially reject events which have a vertex at the 
beam- wall. 

Figures [6?2l and [6731 show the distributions of these two quantities for a sample 
of BB data (i.e. off-peak subtracted from on-peak) from Run 6 overlaid with BB 
MC. The data sample contains the BB events from a sample of 2.81 million on- 
peak events which pass all isBCMultiHadron requirements apart from those on 
PrimVtxdr and PrimVtxdz. The cuts used in hadronic selection (PrimVtxdr < 
0.5 cm, |PrimVtxdz| < 6.0 cm) are far from the sharp peaks of these distributions 
and movements around these cut values have an almost negligible effect on the 
MC time-variation or BB efficiency. 

To ensure the selected events are hadronic, the requirements of BGFMul- 
tiHadron and having at least three good tracks are necessary. There are two 
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E W 




PrimVtxdr (cm) 



Figure 6.2: PrimVtxdr for a data sample from Run 6 (solid line), overlaid with 
BB MC. 




PrimVtxdz (cm) 



Figure 6.3: PrimVtxdz for a data sample from Run 6 (solid line), overlaid with 
BB MC. 
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remaining cut quantities (in addition to plMag) which can be optimized: R2AU 
and ETotal. The total number of cutsets trialled during the optimization pro- 
cedure is 36 X 26 X 12 = 11232. 



Table 6.2: Cut optimization of isBCMultiHadron. 



Quantity 


No. of increments 


Range 


R2A11 


36 


0.45 - 0.80 


ETotal/eCM 


26 


((2.5 - 5.0)/10.58) 


plMag (> 1 DCH hit) 


12 


3.0 - 5.0 GeV/c, or no cut 



The optimization procedure is very similar to the mu-pair optimization. 
There are 66 months of uds and cc on-peak MC, 27 of uds off-peak MC and 
26 of cc off-peak MC. The statistical uncertainties are again treated binomially 
and are given by (|5.3p . The weighted mean efficiency and statistic are given 
by ()5.4|) and ()5.5|) respectively. 

The time variations of these four types of MC are shown in Fig. 16. 41 - 16. 71 The 
cutset pictured has fairly typical time-variation: R2A11 < 0.50, (ETotal/eCM) 
> 0.42533 and plMag < 3.8 GeV/c. These plots all display obvious systematic 
time- variation. The shapes of the distributions are reasonably similar across all 
four plots and they also resemble the muon-pair efficiency variation of Fig. 15.31 
and EH 

The distributions of for all 11232 cutsets are shown in Fig. 16.81 and 16.91 
These are all clearly well above the M — 1 values for each type of MC. 

A comparison between the new cutsets and the existing B Counting selector, 
isBCMultiHadron can be made. The on-peak (off-peak) statistic for this 
tag is 1081.82 (527.27) for uds MC and 2669.32 (1304.13) for cc MC. Almost 
all cutsets for each MC type have values below these which indicates an 
improvement on average. 

6.5 Estimating Systematic Uncertainty 

To estimate the systematic uncertainties in continuum MC time-variation, we 
use the procedure outlined in Section 15.51 We look for in the modified 
statistic of (j5.6p such that 

x"(do) = A/-l. (6.2) 

Again, do is found numerically by incrementing d from zero in steps of 10^^ 
until x'^(d) < M ~ I. The ranges of do values for different cutsets are (1.07 - 
1.69) X 10"^ in on-peak uds MC, (0.67 - 1.59) x lO'^ in off-peak uds MC, (1.60 
- 2.13) X 10"^ in on-peak cc MC and (1.29 - 2.16) x lO'^ in off-peak cc MC. 

As with muon-pairs, the total uncertainty for a particular cutset is then the 
sum in quadrature of statistical uncertainty (l/-y/iy) and systematic uncertainty 
(do)- The 'optimal' set of hadronic cuts is determined at the same time as the 
'optimal' set of muon-pair cuts by finding the combination which causes the 
lowest uncertainty in the total number of B mesons. 
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Figure 6.4: Time-variation of on-peak uds MC efficiency. Details of cutset 
choice arc in text. The vahic is 372.46 and M (the num- 
ber of degrees of freedom) is 66. The weighted mean efficiency is 
0.52376 ± 0.00008. The dashed fines show ±0.5% from the mean. 
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Figure 6.5: Time-variation of off-peak uds MC efficiency. Details of cutset 
choice are in text. The value is 160.71 and M (the num- 
ber of degrees of freedom) is 27. The weighted mean efficiency is 
0.52557 ± 0.00010. The dashed lines show ±0.5% from the mean. 
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Figure 6.6: Time-variation of on-peak cc MC efficiency. Details of cutset choice 

arc in text. The value is 842.20 and M (the number of degrees of 
freedom) is 66. The weighted mean efficiency is 0.62686 ± 0.00007. 
The dashed lines show ±0.5% from the mean. 
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Figure 6.7: Time-variation of off-peak cc MC efficiency. Details of cutset choice 
are in text. The value is 319.67 and M (the number of degrees of 
freedom) is 26. The weighted mean efficiency is 0.62941 ± 0.000010. 
The dashed lines show ±0.5% from the mean. 
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6.6 Variation in k,x 

Using the efSciencies and cross-sections found above, we can evaluate the con- 
tinuum part of K, defined in (j4.8p . The cross-sections of e+e~ — s- uu/dd/ss/cc 
all scale inversely with s (the CM energy squared): 

£xcrx _ 10.542 
'^^ ^ ^ " 10:58^ " 

Note that in this approximation, the product of kx and (defined in (j5.9p ) is 
exactly unity. 

There are other processes in data which can pass the hadronic selection 
and can affect kx but are not included in this approximation. These include 
two photon events. Initial State Radiation, cosmic rays, Bhabhas and beam gas 
events. The hadronic selection is designed to minimize many of these (especially 
cosmic rays and beam gas events) and the effects of these sorts of contamination 
are discussed in Chapter [Sj 

To ensure there is on- and off-peak uds and cc MC in every time bin, the 
time- variation of kx is divided into 23 unequal time periods. As with n^, the 
statistical uncertainty of kx in each time bin includes the statistical uncertain- 
ties of each type of MC given by (|5.3p (e.g. A{euds,i)) and a contribution from 
the uncertainty in the MC cross-sections (e.g. A{(7uds))- 

To an excellent approximation, 

^uds^uds ~^ ^CC^CC { r. \\ 

Kx ; \ — ; — r (o-4) 

uds uds cc cc 

and the statistical uncertainty in kx for each time period is given by (|X6|) : 



A{kx.^) 



e' , g' g'- 

, uas.t uds cc.i cc 



A(ecc,iO-cc) 



uds.i uds 



^i^'uds,i'^uds) ' {Euds,i<^uds + Ecc,i<^c 



i^uds,i^uds + ^cZA^cZ))"^ 



ds 



^ uds,i uds cc,i cc) ' 



where, following (|5.10p . 

A{eudsA<^uds) — [[A.{euds,i)^udsY + {£udsAA{<7uds)T] 



(6.5) 



(6.6) 



and similarly for other MC types. 

The time- variation of kx for the same cutset shown in Fig. 16. 41 - 16. 71 is shown 
in Fig. 16.101 The variation appears to be largely (though not entirely) statistical 
and systematic effects are not as obviously apparent as they are in the individual 
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MC efficiencies. Tlie statistic of kx can be calculated, and for every cutset 
trialled it is above M — 1, so there is a small amount of variation which cannot be 
explained by statistical uncertainties alone. For the cutset pictured in Fig. 16.101 
the statistic is 93.82, and the weighted mean value of kx is 

= 0.98909 ±0.00019, (6.7) 

where the uncertainty quoted is statistical only. The distribution of values 
for all cutsets is shown in Fig. 16.111 Across all cutsets, the mean value of 
is 86.53. The weighted mean values of kx are in the range (0.98751 - 0.99238) 
with a typical statistical uncertainty of 0.00020 or less. 

To estimate the systematic uncertainty in each cutset's kx, again we find do 
satisfying 

x''(do)=M-l, (6.8) 



where 

^^'(^)-E rAr" ^^2l'^2 - (6-9) 
^ {A{Kx,t)r + d^ 

In this case, M - 1 = 22. The do values found are between 0.00066 and 0.00144 
and the distribution of these values is shown in Fig. 16.121 The systematic 
uncertainty in each cutset is added in quadrature with the statistical uncertainty 
to give the total uncertainty on each cutset's kx value. 
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Figure 6.8: Distribution of statistic for trialled cutsets (left) for uds on-peak 

MC and (right) for uds off-peak MC. The mean vahies of these 
distributions are 479.86 and 232.47 respectively. M is 66 and 27 
respectively. 
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Figure 6.9: Distribution of statistic for trialled cutsets (left) for cc on-peak 
MC and (right) for cc off-peak (right) MC. The mean values of these 
distributions are 1292.82 and 441.28 respectively. M is 66 and 26 
respectively. 
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Figure 6.10: Time- variation of kx for the cutset shown in Fig. 16.41 - 16.71 The 

value is 93.82, M is 23 and the weighted average is 0.98909. The 
dashed lines show ±0.5% from the mean. 
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Figure 6.11: Distribution of kx statistic for all trialled cutsets. M is 23. 
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Figure 6.12: Distribution of Kx do values for all trialled cutsets. 
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Chapter 7 

Cutset Optimization 



The main purpose of the research described in this thesis is to recommend an 
optimal set of cuts which count B meson events with a minimal systematic 
uncertainty. This uncertainty has contributions from the uncertainties in the 
counting of both muon-pairs and hadronic events. The number of B meson 
events counted with the optimal cutsets should also agree with the original B 
Counting code, which has been in place since 2000 and has a total uncertainty 
of ±1.1%. 

As was mentioned in Chapter^ the B Counting uncertainty for optimization 
purposes found in this chapter is different to the final estimate of systematic 
uncertainty in data (given in Chapter [S]). In this chapter we find the B Counting 
cutsets which give the lowest sensitivity to the time- variation of wds, cc and mu- 
pair Monte Carlo efficiencies. We also include the statistical uncertainty of the 
MC cross-sections (included in k uncertainties), 

The method used to calculate the propagation of uncertainties is outlined in 
Appendix [Xj 

7.1 B Counting Uncertainty 

The B Counting formula of (|4.17p is: 

K = ^{Nh-N^-RoB-k) (7.1) 
^ {Nh - Km • «x • Rt. ■ N'h) (7.2) 



where 



^ (7.3) 



N' 



N'jj is the number of off-peak hadronic events counted and and kx are defined 
in (|4.8[) . The subscripts "i?" and "X" in an off-peak sample are interchangeable 
since the off-peak hadronic data contains no BB events from T(4S') decays. 
The uncertainty in the number of B meson events, A(A'^^) is given by: 



(7.4) 



Here, Xi represents each quantity with associated uncertainty (such as Suds or 
K^). To estimate the uncertainty in the number of B meson events in a sample of 
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on-peak data of luminosity C with continuum subtracted by an off-peak sample 
of luminosity we can rewrite (|7.f P as: 

K = —[{H + Ho)-K^KxR,AH' + H'o)] (7.5) 



SB 



where 
and 



H = C{euds(^uds + eccCTcc + sbc^b) (7.6) 

Ho and (the subscript "O" represents "other" event-types) are the numbers 
of on- and off-peak events which are not BB^ uds or cc but pass the hadronic 
cuts. These terms are discussed in detail in Chapter HI but are neglected during 
this optimization procedure. We can also write 



where Mq (Mq) is the number of on-peak (off-peak) events which are not 
e~^e~ IJ-^IJ-~ but still pass the muon-pair cuts. We again neglect these during 
the optimization. 

In calculating the uncertainty in the number of B mesons, we omit the 
uncertainties in efficiency of off-peak MC. The time- variations of on- and off- 
peak MC are correlated, and it is important not to double-count the uncertainty. 
Off-peak data (and hence also MC) is taken only for short periods of time around 
longer periods of on-peak running. For this reason, off-peak running can be 
viewed as a 'snapshot' of the non-BB continuum at a particular time. 

It is clear the two are not independent, since and kx both have 
statistics much lower than any MC type has alone. The time-variation visible 
in off-peak MC is to a large degree the same as that displayed in on-peak MC. 
For these reasons, the off-peak MC efficiencies are treated as constants in (|7.4|) . 

Hence, we have: 



' ""-^Kx + tt^Ak^ (7.9) 



where 



\ Okx J V 



{^Eudsf = (A(e„d.)stat)^ + ida.udsf (7.10) 



and similarly for e^c, ^^"^ '^m- Equation (|7.f Op gives the combined sta- 

tistical and systematic uncertainty found during the previous optimization pro- 
cedure and has the same form as the denominator of the x'^ statistic in (|5.6p . 
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Once an optimal cutset is found, a value for Aeb can be estimated by com- 
paring BB MC to data. During the optimization procedure, we estimate its 
value by setting: 

Aeb = 0.2(1 -£b). (7.11) 

This is chosen to favour cutsets with higher BB efficiency. For a cutset with is 
of 0.96, (|7.1ip gives an uncertainty in is of around 0.83%. 
The differential quantities of (|7.9p are: 



dN^ 1 



dN'^ _ 1 



—Ca^ds (7.12) 



(7.13) 



*7 ^ -rs^'-^-'e±:" 



dux £b 



Ht.Rt.H' (7.15) 



dN% 1 



— KxRi^H' (7.16) 



and 



dNl _ 1 
den 



-2 



(-^ ^£uds^uds £cc^cc) 

KRf,£' {e'^as<^'^ds + ^Wcc) ) • (7-17) 

This last term is very small (and would be exactly zero if B Counting were 
perfect) so is negligible in comparison to the other terms. 

As was mentioned above, the expression (j7.9p is useful for finding an optimal 
set of B Counting cuts, but it will not be used to estimate the systematic 
uncertainty of data. 



7.2 Optimization 

The optimal mu-pair and hadronic B Counting cutsets are those which minimize 
{AN'^)/Ng (the relative uncertainty in Ng) while meeting the selection goals 
outlined in Section [47T] As a cross-check, the number of B mesons counted in a 
sample by any pair of cutsets can be compared to that given by the existing B 
Counting code. 

The total uncertainty of ()7.10|) can be divided into two parts: one depen- 
dent only on uncertainties related to muon-pair selection and the other only on 
hadronic uncertainties. The muon-pair part 
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includes terms dependent on Ae^ and Ak^. The hadronic part A{N'^)'jf con- 
tains all other terms. The two are not entirely independent — for example the 
Ak^ term, (|7.16p depends specifically on the number of off-peak hadronic events 
counted by a chosen set of hadronic cuts. 

The only quantities in (|7.9|) which are constant independent of cutset choice 
are £, £' and the cross-sections. The variable quantities of (|7.12| - [7?TB)) — the 
number of hadronic and muon-pair data events, efficiencies, kx and are 
specific to each cutset trialled during the procedure. 

The first stage of the optimization procedure calculates the number of B 
mesons in a sample of on-peak data and the value of A(A''^)^ for each muon- 
pair cutset. Continuum subtraction is provided by a corresponding off-peak 
data sample. 

To calculate these values, the number of hadronic events in the two samples 
must be input. Initially these are taken to be the number of isBCMultiHadron 
events in each sample, i.e. the number counted by the existing B Counting code. 
It is also necessary to input the value of es for isBCMultiHadron (approximately 
0.96) and a value of kx, which is initially taken to be the theoretical estimate 
of ((O)) . The 61880 muon-pair cutsets are then sorted by A{N%)f,/N^. 

Assuming it meets the B Counting selection goals, the 'optimal' muon-pair 
cutset is the one with the lowest relative uncertainty in Ng. The numbers of 
muon-pair on- and off-peak data events counted by this cutset, along with the 
value of then used as inputs to optimize the hadronic selection. 

Following the same procedure, the values of N'^ and A(A^^)|^ are calculated 
for the 11232 hadronic cutsets, which are then sorted by A{Ng)'jf /Ng. The 
values of Eb, kx and the number of on- and off-peak hadronic events counted 
by the optimal hadronic cutset are used as inputs to re-optimize the muon-pair 
cutsets. This procedure continues until the optimal cutset choice stabilizes. 

This procedure was performed on each major running period (Runs 1-6). 
There is typically less than 5% variation in the value of A{Ng)fj_/Ng across 
the first five hundred muon-pair cutsets (this number is around 10% for the 
same number of hadronic cutsets). With so many nearly-identical cutsets, the 
optimal choice for each Run varies and statistical fluctuations can change each 
cutset's rank. 

Due mainly to these statistical fluctuations, no single cutset is consistently 
the optimal choice, but some trends emerge. Muon-pair selection cutsets tend 
to have a low values of A{Ng) / Ng if: 

• the masspair/eCM cut is high, typically in the range (0.756 - 0.827), 

• the maxpCosTheta cut is high, typically above 0.735, 

• the acolincm cut is high, typically above 0.11. 

No trend is obvious for the maxpEmcLab cut. For the hadronic selection, cutsets 
tend to have low values of A{N^)h/N'^ if: 

• the ETotal/eCM cut is high, typically above 0.378, 
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• the plMag cut is low, typically less than 4 GeV/c. 

No trend is obvious for the R2 cut, and cutsets without a plMag cut do not 
generally perform near the top. 

7.3 Choice of Muon-Pair Selection Cutset 

As mentioned above, many cutsets have almost identical values of A{Ng)fj,/Ng 
or A{Ng)if/Ng. The change to the new B Counting particle lists (Charged- 
TracksBC and NeutralsBC) has increased stability (in comparison to isBC- 
MuMu and isBCMultiHadron) even before any cuts were varied. Deciding on 
one cutset to recommend for B Counting is to some degree an arbitrary choice. 
The cutset we recommend to select muon-pairs is: 

• BGFMuMu, 

• masspair/eCM > 8.25/10.58, 

• acolincm < 0.13, 

• uTracks > 2, 

• maxpCosTheta < 0.745, 

• either maxpEmcCandEnergy > GeV, or at least one of the two highest- 
momentum tracks is identified as a muon in the IFR, 

• maxpEmcLab < 1.1 GeV. 

This can be compared to the isBCMuMu requirements listed in Table [4711 Note 
that the meaning of nTracks is slightly changed, as it now refers to the number 
of tracks on the ChargedTracksBC list. 

Some properties of this cutset and isBCMuMu for muon-pair MC are given 
in Table O 

This new cutset clearly displays less systematic variation than isBCMuMu. 
For example, the difference in c?g values shows an order of magnitude stability 
improvement in k^. The time- variation of on- and off-peak MC for these new 
cuts is shown in Fig. 17.11 and 17.21 The equivalent plot for the isBCMuMu tag 
(with a vertical axis of the same scale) is shown in Fig. l7.3l and l7.4l and the time 
variations of are shown in Fig. 17.51 and 17.61 

We use one other method to estimate the systematic uncertainty on — we 
vary the mu-pair selection cuts on maxpCosTheta and measure how the change 
affects the quantity 

N 

«M • 7^ = ^^R^- (7-19) 

This quantity appears explicitly in the B Counting formula (|4.17p . Any changes 
to it as the selection cuts vary show the combined effects of the altered cuts on 
both data and MC. This allows us to estimate the systematic uncertainty on 
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as the number of mu-pair events counted should be almost constant regardless 
of the arbitrary choice of mu-pair selection cuts. 

Specifically we allow the cuts on maxpCosTheta to vary between 0.60 and 
0.77 (in increments of 0.005) while the value of K^^i?^ is recalculated for each cut 
choice. The value of n^Rfj, for the original cuts is 9.7126. The largest (positive 
or negative) variation from the original value occurs when maxpCosTheta is 
equal to 0.635 — this causes a decrease in K^i?^ of 0.047%. We attribute 
this percentage difference to an additional systematic uncertainty in as an 
estimate. 



Table 7.1: Properties of the recommended muon-pair cutset and isBCMuMu. 





Cutset property 


New Cutset 


isBCMuMu 


On-peak: 




0.41558 


0.42566 




Stat, uncert. in 


0.00010 


0.00010 




statistic of 


118.11 


1016.58 




Syst. uncert. in {do n) 


0.00082 


0.00311 


Off-peak: 




0.41598 


0.42627 




Stat, uncert. in e' 


0.00014 


0.00014 




statistic of 


51.71 


402.25 




Syst. uncert. in e'^ ((io,^) 


0.00060 


0.00353 




^/i 


1.00791 


1.00769 




Stat, uncert. in 


0.00044 


0.00043 




statistic of 


23.82 


62.11 




Syst. uncert. in (rfg,^) 


0.00032 


0.00306 




Syst. uncert. in {k^R^) 


0.00047 


N.A. 



7.4 Choice of Hadronic Selection Cutset 

The cutset we recommend to select hadronic events is: 

• BGFMultiHadron, 

• R2A11 < 0.65, 

• nTracks > 3, 

• ETotal/eCM > 4.60/10.58, 

• PrimVtxdr < 0.5 cm, 

• |PrimVtxdz| < 6.0 cm, 

• plMag < 3.0 GeV/c (> 1 DCH hit). 
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Figure 7.1: On-peak muon-pair MC efficiency variation for the new proposed 
muon-selection cuts. The dashed lines show ±0.5% from the mean. 
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Figure 7.2: Off-peak muon-pair MC efficiency variation for the new proposed 
muon-selection cuts. The dashed lines show ±0.5% from the mean. 
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Figure 7.3: On-peak muon-pair MC efBciency variation for isBCMuMu. The 
dashed hnes show ±0.5% from the mean. 
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Figure 7.4: Off-peak muon-pair MC efficiency variation for isBCMuMu. The 
dashed lines show ±0.5% from the mean. 
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Figure 7.5: Time-variation of for the new cutset. The dashed hnes show 
±0.5% from the mean. 
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Figure 7.6: Time-variation of for isBCMuMu. The dashed Unes show ±0.5% 
from the mean. 
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Some properties of this cutset and isBCMultiHadron are shown in Table 
17.21 The proposed cutset is more stable in time than isBCMultiHadron, but 
the improvement from the existing tag is not as large as it was for the muon- 
pair selector. Nevertheless the new hadronic cutset has a lower statistic for 
all types of MC, and consequently the estimated systematic uncertainty (do) is 
lower in each case. The statistics of each MC type decrease to between 26% 
and 36% of their isBCMultiHadron values. 

The time-variation of nx and the efficiency of each type of continuum MC 
for the proposed cutset and isBCMultiHadron are shown in Fig. 17.91 - 17.181 The 
time- variation of uds MC efficiency for the new cuts and isBCMultiHadron are 
shown in Fig. 17. 91 - 17. 121 The equivalent plots for cc MC efficiency (with vertical 
axes of the same scale) are shown in Fig. 17.131 - 17.161 Finally, kx for the new 
cutset and isBCMultiHadron is shown in Fig. 17.171 - 17.181 

The weighted mean BB efficiency (es) for this cutset is 0.94046 ± 0.00003 
(where the uncertainty is statistical only), which is slightly lower than the is- 
BCMultiHadron value of 0.9613. The main reasons for this decrease are the 
addition of a plMag cut and the slightly higher cut on ETotal. 



7.4.1 ISR Production and Two-Photon Events 

We have not yet considered the effect on kx from events which are difficult to 
simulate with Monte Carlo; in particular Initial State Radiation (ISR) and two- 
photon events. Our limited knowledge of these events introduces an additional 
systematic uncertainty to the value of kx- 

ISR events are those in which one of the incoming beam particles emits a 
photon before interacting. The theoretical ISR cross-section used in this analysis 
is based on that described by Benayoun et al. [25]. For ISR production of any 
state (labelled "V"), the ISR cross-section is taken to be: 

where is a function defined in [25j . and my and T^e are the resonance mass 
and partial width oiV ^ e+e~ respectively. The total cross-section is estimated 
to be the sum of each possible J^'~^ = 1~~ resonance. The dominant terms in 
this sum are to the p° (by an order-of-magnitude), J/iji and T(3S') resonances. 

Two-photon events are those where the two colliding particles emit virtual 
photons which in turn interact to produce some system of particles. A review 
of two-photon particle production can be read in [26j . For high energies (such 
as those provided by PEP-II), the cross-section for two-photon production of a 
lepton pair is given approximately by: 

28a4 / s \\ 2 

ae+e-^e+e-^+£- = ^ ] ^ ^ (7-21) 



where / can be either e or /i and ^/s is the CM energy. 
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The (high-energy hmit) cross-section for hadron production (h) is given ap- 
proximately by: 

sml sm^ f s 

(Te+e-->e+e-h= .r, 22 '^^ 2 2 ^" 6 2 ^ ' C^'^^) 

IHtt^to^ m^mi m^m^ \ m%. J 

where Wp and m-^ are the masses of the p and tt mesons. 

Since these two processes are not known especiahy weh, we attempt to quan- 
tify their effect on systematic uncertainty. We do this by using cc, uds, t^t^ 
and high-angle Bhabha MC to estimate the effective cross-section for other (i.e. 
not of these four types) events in off-peak data of passing the hadronic selection, 
and the corresponding systematic uncertainty on k. 

The hadronic part of k takes the form: 

KX = -r—r (7-23) 



E known , \-^unknown 

E known t i , \-^unknown i > 



(7.24) 



Here, by known we mean events for which BaBAR has a reliable MC model 
and well-known cross-sections. The effects of the unknown types can be esti- 
mated by comparing off-peak data and the known MC. We assume that the un- 
known component is comprised solely of two-photon, ISR and low-angle Bhabha 
events. 

We compare the known MC types to off-peak data. The cross-sections for 
each MC type are used to combine them in the correct proportions to model 
data as accurately as possible. Normalization is achieved by forcing the MC to 
have the same value as data at 0.9 on the plot of total event energy. The MC 
distributions of total energy are shown in the left plot of Fig. 17.71 and the sum 
overlaid with data on the right. 

The excess in data over the cut range (above 0.4348 on the ETotal/ y^{s) 
plot) is found to be equivalent to an effective cross-section of approximately 
0.10 nb. The value of kx for the optimal cuts for only the known types of 
MC has the value 0.98979 ± 0.00016 (where the uncertainty is statistical only). 
To estimate kx including two-photon and ISR events we estimate the effec- 
tive on-peak cross-sections by scaling by the theoretical cross-sections. These 
considerations increase the value of kx and introduce an additional systematic 
uncertainty on its value. We combine this new uncertainty in quadrature with 
that due to MC time-variation. 

Note that we do not use the (fairly old) theoretical models to predict the 
absolute cross-section, but use the ratio of theoretical on- and off-peak cross- 
sections for the processes to estimate the total effective on-peak cross-section. 

The proportions of two-photon and ISR events in the data are unknown, so 
for a range of a values from 0.01 to 0.99 we let 

unknown 

^ e'jaj = a ■ e'jsj^a'jsj^ + (1 - q) • £27^^27- (7-25) 
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and use the ratios of (^isrI '^'isb. ^'^'^ '^'i-i I '^'i-i from theory to estimate the on- 
peak unknown effective cross-section (assuming the on- and off-peak efficiencies 
are identical). 

The value of kx including the unknown component can be estimated. We 
take its value to be the a = 0.5 value, and equate the difference as a varies to 
be systematic uncertainty. Performing the identical analysis for the theoretical 
Bhabha cross-section does not change the solution — for 100% Bhabhas or a 
1:1:1 mixture of two-photon, ISR and Bhabha events, the value of kx stays 
within the range found using two-photon and ISR events alone. 
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Figure 7.7: Comparison of ETotal/yi for known MC types to off-peak data for 
events passing all other hadronic selection cuts: (left) in decreasing 
order distributions of sum of all MC types, uds, cc, T"''r~, low-angle 
Bhabha MC; (right) the sum (points) is shown overlaid with data - 
the plots are normalized to have the same value at 0.9. 

Fig. 17.81 shows the ETotal/eCM distribution of the "unknown" events (i.e. 
the difference between the MC and off-peak data distributions). 

The resuh is: kx = 0.9901 ± 0.0002. Note that the addition of two-photon 
and ISR events give a central value of kx higher than that calculated for the 
known types only (although they agree within uncertainties). We assume this 
uncertainty is independent of that calculated in Section [6.61 

Combining the two uncertainties (from the MC time- variation and this two- 
photon, ISR and low-angle Bhabha analysis) in quadrature, the value for the 
optimal cutset is: kx = 0.9901 ± 0.0010. 



Chapter 7. Cutset Optimization 65 



90000 
80000 
70000 
60000 
50000 
40000 
30000 
20000 
10000 


0.2 0.4 0.6 0.8 1 1.2 

Total Event Energy/sqrt(s) 



Figure 7.8: Difference in ETotal between known MC types and off-peak data. 
Tfie cut at 0.4348 is marked with a dashed line. 



7.5 B Counting Results 

The numbers of B meson events counted in data by this pair of new proposed 
cutsets for each major running period are given in Table 17.31 The two uncer- 
tainty terms, A{Ng)f^/N^ and A{N^)h /N'^ are added in quadrature to give 
the total uncertainty in the number of B mesons. For comparison, the number 
counted by the existing code (with a total uncertainty of 1.1% in each case) is 
also provided. 

For every Run, the number of B meson events counted by the new cutsets 
agrees with the existing value within uncertainty. The total uncertainty ranges 
between 0.615% and 0.634%. Note that these values are calculated only to find 
the optimal pair of cutsets. The final estimate of systematic uncertainty in data 
using the optimal cutset is discussed in Chapter [51 
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Figure 7.9: On-peak uds MC efficiency time-variation for the new proposed 
hadronic cuts. The dashed hnes show ±0.5% from the mean. 
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Figure 7.10: Off-peak uds MC efficiency time-variation for the new proposed 
hadronic cuts. The dashed lines show ±0.5% from the mean. 
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Figure 7.11: On-peak uds MC efRciency time- variation for isBCMultiHadron. 
The dashed lines show ±0.5% from the mean. 
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Figure 7.12: Off-peak uds MC efficiency time-variation for isBCMultiHadron. 
The dashed lines show ±0.5% from the mean. 



Chapter 7. Cutset Optimization 



68 



0.762 



E 0.76 



0.758 
0.756 
0.754 
0.752 
0.75 
0.748 



>_ > > I > 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 






'l''!'" ,||",V ^ 








: |t : 

7 , , 1 , , , 1 , , , 1 , , , 1 , , , 1 , - 



20 



40 



60 



80 100 
Month number 



Figure 7.13: On-peak cc MC efficiency time- variation for tfie new proposed 
hadronic cuts. The dashed hnes show ±0.5% from the mean. 
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Figure 7.14: Off-peak cc MC efficiency time- variation for the new proposed 
hadronic cuts. The dashed hnes show ±0.5% from the mean. 
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Figure 7.15: On-peak cc MC efficiency time- variation for isBCMultiHadron. 
The dashed lines show ±0.5% from the mean. 
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Figure 7.16: Off-peak cc MC efficiency time-variation for isBCMultiHadron. 
The dashed lines show ±0.5% from the mean. 



Chapter 7. Cutset Optimization 



70 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
10 12 14 16 18 20 22 
Time bin number 



Figure 7.17: Time- variation of kx for the new cutset. The dashed lines show 
±0.5% from the mean. 
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Figure 7.18: Time- variation of kjc for isBCMultiHadron. The dashed lines show 
±0.5% from the mean. 
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Table 7.2: Properties of the recommended hadronic cutset and isBCMulti- 
Hadron. 





J. J. 1. 

Cutset property 


New Cutset 


isaC- 

1\ /Till + i U.i -^-..^-.1-1 

IVlmtiJrLaQi on 


On-peak: 


^uds 


u.D/yys 


U.Dl / Do 




St3-t. uiiCGrt. in Suds 


u.uuuu / 


u.uuuu / 




statistic of Suds 


344.93 


1081.82 




Syst. uncert. m iuds {dQ^uds) 


0.00121 


0.00227 






0.75364 


0.73431 




Stat, uncert. in Scc 


0.00006 


0.00006 




statistic of e^- 


939.95 


2669.32 




Syst. uncert. in e^c (do,cc) 


0.00175 


0.00278 


Off-peak: 


p' 

^uds 


0.63200 


0.61752 




Stat, uncert. in e^^^ 


0.00009 


0.00009 




statistic of e'^j^^ 


164.46 


527.27 




Syst. uncert. in i'^^^ (rfo,«dJ 


0.00094 


0.00216 






0.75568 


0.73580 




Stat, uncert. in e'^- 


0.00009 


0.00009 




statistic of e^^ 


349.56 


1304.13 




Syst. uncert. in e'^- (cIq^^-^) 


0.00167 


0.00274 




kx {known MC only) 


0.98979 


0.99198 




Stat, uncert. in kx 


0.00016 


0.00017 




statistic of kx 


80.19 


101.19 




Syst. unccrt. in k y {d()roiii) 


0.00095 


0.00122 




S>'sl. unccrt. in (ISR./2";) 


0.0002 


0.0002 




kx (incl. ISR/27) 


0.9901 


0.9918 




Combined uncert. in kx 


0.0010 


0.0012 
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Table 7.3: Uncertainty terms and number of B meson events with overall un- 
certainty found during optimization procedure in Runs 1 6. 





Run 1 


Run 2 


N'^ (new cutsets) 


0.461% 

0.416% 
(22.39 ±0.14) X 10^ 


0.456% 
0.412% 
(67.61 ±0.42) X lO'* 


(existing code) 


(22.40 ±0.26) X lO'* 


(67.39 ±0.74) X lO*" 




Run 3 


Run 4 


(new cuts(>tri) 


0.463% 
0.417%) 
(SS.?.-) ± 0.22) X lO" 


0.460% 
0.414% 
(110.70 ± ().(i9) X 10'' 


Ng (existing code) 


(35.57 ±0.39) X 10'* 


(110.45 ± 1.21) X 10" 




Run 5 


Run 6 


Ng (new cutsets) 


0.462% 
0.416% 
(146.19 ±0.91) X 10^ 


0.471% 

0.424% 
(82.19 ±0.52) X 10'* 


A^)^ (existing code) 


(147.19 ± 1.62) X 10** 


(82.04 ±0.90) X 10** 
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Chapter 8 

Sources of Contamination 



As mentioned in previous Chapters, there are several possible contamination 
effects which are considered negligible for B Counting. This includes e~^e~ 
T~^T~ decays incorrectly passing hadronic selection and Bhabha events incor- 
rectly passing muon-pair selection among others. In this Chapter we quantify 
these effects. 



8.1 Bhabha Events 

The process with by far the highest cross-section in PEP-II is Bhabha scattering. 
Most processes, like e+e~ M^M" or e+e" — > qq have cross-sections propor- 
tional to where Ecm is the event's CM energy. The cross-section for 
Bhabha scattering is more complicated, and at PEP-II's on-peak CM energy of 
10.58 GeV is around 40 times larger than most other processes. The differential 
cross-section takes the form: 



da 
d cos 



1 1\^ /tV fs"-^ 



\ - + - 



t \s \t 



(8.1) 



where s, t and u are the Mandelstam variables, which relate the four-momenta 
of the incoming and outgoing particles [22] . 

The majority of Bhabha events scatter at low angles (i.e. continue down the 
beam-pipe), but the cross-section for large-angle Bhabha scattering (i.e. within 
the angular limits of ChargedTracksBC) is still substantial. Based on samples 
of one million generated Monte Carlo Bhabha scattering events, the values for 
on- and off-peak running are (13.8910 ± 0.0040) nb and (14.0918 ± 0.0041) nb 
respectively. These are more than ten times the BB cross-section of 1.1 nb. 

Ideally, the B Counting muon-pair and hadronic selectors should reject all 
Bhabha events, but due to many imperfections a small number incorrectly pass. 
To estimate the size of this effect, we apply the selectors to 7.725 million high- 
angle Bhabha MC events. The numbers of Bhabha MC events passing the 
selectors are given in Table lOI For the new cuts, the Bhabha hadronic fake-rate 
is approximately 0.006%, and the Bhabha muon-pair fake-rate is approximately 
0.002%. 

We can estimate the ratio of Bhabha events to B meson events passing the 
hadronic cuts in an on-peak sample of luminosity C: 



Nsha (pass had, cuts) _ esha'^BhaC 
Nb (pass had. cuts) e^'^ asC 



(8.2) 
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Table 8.1: Fake-rates for Bhabha MC events passing hadronic and mu-pair 
selectors. The sample contained 7.725 million large-angle on-peak 
Bhabha MC events. 





Number passing selector 


Percent of total 


New hadronic selector 


449 


0.0058% 


isBCMultiHadron 


1559 


0.0202% 


New mu-pair selector 


149 


0.0019% 


isBCMuMu 


22 


0.0003% 



Nsha (pass had. cuts) _ 0.00006 x 13.9 nb 
Nb (pass had. cuts) ~ 0.941 x 1.1 nb ^ ' ^ 

w 0.08%. (8.4) 

Thus, even though the cross-section for high-angle Bhabha events is more than 
ten times the other continuum cross-sections, the Bhabha hadronic fake-rate is 
very much smaller than the selector's systematic uncertainty. The equivalent 
number for the muon selector (of fake Bhabhas to genuine muon-pairs) is factor 
of three smaller. 

The final possible source of Bhabha contamination relates to the energy- 
dependence of the Bhabha cross-section. One of the most important assump- 
tions for B Counting is that background events scale in the same way with CM 
energy as muon-pair events. This assumption enables backgrounds to be sub- 
tracted by using the ratio of on-peak to off-peak muon-pair events. For Bhabha 
events however, this assumption is not true. The ratio of on-peak to off-peak 
muon-pair MC events from Section ISTTl is: 



and for Bhabhas, this ratio is 

OBha ^ 13.8910 nb 
14.0918 nb 



= 0.98575 (5 sig. fig.). (8.6) 



This difference means that when subtracting continuum events, the number 
of on-peak Bhabha events passing the hadronic cuts will be over-estimated by a 
factor of approximately 0.7%. The total effect is 0.7% of the 0.006% of Bhabha 
events passing hadronic selection. Since this number is negligibly small it can 
be ignored safely. 



8.2 Continuum Fake-Rates 



There is a small probability of hadronic continuum events passing muon-pair 
selection. These events scale with energy in the same way as muon-pairs. When 
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counting B meson events, the ratio of muon-pairs (i?^) is used, so if the per- 
centage of on- and off-peak fake muons from hadronic continuum events are the 
same, any hadronic continuum fakes will have no effect on B Counting. 

Table [8T2l shows the percentage of hadronic continuum MC events which pass 
the new muon-pair selector and isBCMuMu. Clearly uds events have a higher 
fake-rate, but the overall effect is too small (in both the number of fakes and 
on-/off-peak difference) to have any consequences for B Counting. 

Table 8.2: Fake-rates for hadronic continuum MC events passing the mu-pair 



selector. The sample contained 5.20 million on-peak cc, 4.20 million 
off-peak cc, 6.40 million on-peak uds and 4.20 million off-peak uds 
MC events. 





% passing new muon-pair selector 


% passing isBCMuMu 


cc (on-peak) 
cc (off-peak) 


0.00014% 
0.00002% 


0.00036% 
0.00021% 


uds (on-peak) 
uds (off-peak) 


0.009% 
0.010% 


0.024% 
0.026% 



There is also a small chance that muon-pair events can pass the hadronic 
selector. The probabilities for this are shown in Table [8T3l The effect is negligibly 
small. 



Table 8.3: Fake-rates for muon-pair MC events passing the hadronic selector. 

The sample contained 4.05 million on-peak and 4.08 million off-peak 
/i+Zi" MC events. 





% passing new 
hadronic selector 


% passing isBCMulti- 
Hadron 


(on-peak) 
(off-peak) 


0.0078% 
0.0082% 


0.0359% 
0.0377% 



8.3 Tau Production 

Another process with a cross-section close to 1 nb is e+e^ t^t^. The r 
leptons are short-lived and decay within the beam pipe. These events have only 
a small probability of passing the muon-pair selector, but the probability of 
passing the hadronic selector is much larger. 

Similarly to mu-pair events, tau-pair events are simulated with the KK2F 
MC generator, GEANT 4 and background addition. The production cross- 
sections for the process e+e^ t^t~ based on 10 million generated on- and 
off-peak events are (0.918797 ± 0.000066) nb and (0.925180 ± 0.000066) nb re- 
spectively. 

The hadronic and mu-pair selector efhciencies are shown in Table [831 The 
hadronic efficiency for r pairs is approximately 4% of the uds and cc hadronic 
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efficiencies. 

Table 8.4: Fake-rates for r+r^ MC events passing hadronic and mu-pair se- 
lectors. The sample contained 4.55 million on-peak and 4.26 million 



off-peak e+e 


T^T^ MC events. 






% passing new 


% passing isBCMulti- 




hadronic selector 


Hadron 


r+r^ (on-peak) 


2.67% 


1.37% 


T+T^ (off-peak) 


2.68% 


1.37% 




% passing new 


% passing isBCMuMu 




muon-pair selector 




T^T^ (on-peak) 


0.052% 


0.126% 


r+T^ (off-peak) 


0.05f% 


0.117% 



The rate at which r leptons pass the hadronic cuts is a factor of two higher 
than for isBCMultiHadron, but this largely does not affect B Counting. The 
increase in efficiency is mainly due to the increased R2 cut (which changed from 
0.5 to 0.65). The cross-section for r production changes with energy identi- 
cally to muon-pairs, and since the on-peak and off-peak efficiencies are almost 
identical, the overall effect on B Counting is small. 

As for other MC types, we can calculate the statistic for the time- variation 
of r+T" MC. Again, we find the do value of (|5.6p such that 

x''(do) = M-l. (8.7) 

This value provides an estimate of the systematic uncertainty in the r lepton 
efficiency. 

The time- variation of on- and off-peak r+r^ MC is shown are Fig. 18. II and 
18.21 The overall uncertainties (combining systematic and statistical uncertain- 
ties in quadrature) are ±1.5% for on-peak and ±0.99% for off-peak t lepton 
MC. When combined with the efficiency for passing the hadronic cuts of around 
2.7%, the contribution to B Counting uncertainty from variations in r lepton 
efficiency is an order of magnitude lower than the uds and cc contribution. 
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Figure 8.1: On-peak tau lepton MC efSciency variation for the new proposed 
hadronic cuts. The dashed lines show ±0.5% from the mean. 
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Figure 8.2: Off-peak tau lepton MC efficiency variation for the new proposed 
hadronic cuts. The dashed lines show ±0.5% from the mean. 
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Chapter 9 



Systematic Uncertainty 



In this chapter we estimate the total systematic uncertainty on B Counting. 
We consider the uncertainties on BB efficiency, kx and and estimate the 
combined uncertainty. 

9.1 BGFMultiHadron and Tracking Efficiency 

One contributor to systematic uncertainty in is the efficiency of the BGFMul- 
tiHadron tag. To examine this, we studied samples of random on- and off-peak 
events which pass the Level Three trigger (the stage at which BGF tags are 
assigned). This dataset has the name DigiFLSOpen and consists of 0.5% of 
all Level Three events. The samples contained 17.57 million on-peak and 1.67 
million off-peak data events from Run 4. 

BGFMultiHadron requires at least three charged tracks and R2 of less than 
0.98 in the CM frame. For events with a low number of charged tracks, the 
quantities BGFMultiHadron and nTracks are closely correlated. Events which 
fail hadronic selection because nTracks is too low (less than three) also usually 
fail because they do not have the BGFMultiHadron tag. In the on-peak sample, 
1.89 million events passed all hadronic cuts except those on BGFMultiHadron 
and nTracks. The combination of BGFMultiHadron and nTracks cuts rejected 
a further 0.19 million events; of these, 86% failed both cuts. This number was 
87% in the off-peak data. 

To find the component of systematic uncertainty due to the combination of 
the BGFMultiHadron and nTracks cuts in hadronic selection, we compare the 
low-end tail of the nTracks distribution in BB data and MC for events passing 
all hadronic cuts except those two. The distribution for data (luminosity-scaled 
off-peak subtracted from on-peak) overlaid with BB MC is shown in Fig. 19.11 
The BB MC sample contained 12.43 million events. 

After scaling, the off-peak data has more zero- and one-track events than the 
on-peak data. This causes the equivalent data bins in Fig. 19.11 to be negative. 
This is due in a large part to contamination by Bhabha events. The proportion 
of on-peak (off-peak) data events in this plot with exactly one charged track, 
where that track is identified as an electron is 83.7% (83.8%). This compares 
to 12.1% for BB MC. Similarly the proportion of two-track events where both 
are identified as electrons is significantly higher in data (26.2% and 28.3% in 
on-peak and off-peak respectively) than in MC (3.6%). 

Events with one or two charged tracks where the number of charged tracks 
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Figure 9.1: Distribution of nTracks for BB data (solid) and MC (points). All 
hadronic cuts except those on BGFMultiHadron and nTracks are 
applied. The two are normalized to have the same area from 3-20. 



is equal to the number of electrons in the event are likely to be Bhabhas. When 
tracks of this type are removed, the agreement between data and MC improves, 
especially in the bin of data events with exactly one track. The revised distri- 
bution is shown in Fig. 19.21 The MC distribution is largely unchanged. 

To estimate the contribution to systematic uncertainty by these two cuts 
we can compare the proportion of events in Fig. 19.21 with nTracks less than 
three for BB data and MC. For data we sum the absolute value of the bin 
contents (since they are negative for the first two bins) and find 0.56% of all 
events in the plot are in the first three bins. For BB MC, this number is 0.20%. 
The difference (0.36%) is attributed to systematic uncertainty in eb from the 
BGFMultiHadron and nTracks cuts. 



9.2 Comparison of Data and MC 

Next we estimate the systematic uncertainty in eb for the optimal cutset choice 
due to the difference between BB MC and data. We assumed in Chapter 2] that 
the proportions of charged and neutral BB decays of the T(45') are the same, 
i.e. that 

^^fo _ BiTjiS) ^ B^B_-) _ ^ 

B{r{iS) B"B^) ■ ^ ' ' 

Since the decays of charged B mesons are slightly different to those of neutral 
ones, a source of uncertainty arises if this assumption is false. 
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Figure 9.2: Distribution of nTracks for BB data (solid) and MC (points) with 
'obvious' Bhabha events removed. All hadronic cuts except those 
on BGFMultiHadron and nTracks are applied. In addition events 
are not plotted if the number of electrons is non-zero and equal to 
the number of charged tracks. The two are normalized to have the 
same area from 3-20. 



Two of the most precise measurements of foo = B{T{4:S) B^B°) and 

/_! = B{T{4:S) B'^B~) are described in [27j. That analysis uses a partial 

reconstruction of the decay B'~' D*'^t~v to obtain: 

/oo = 0.487 ± O.OlO(stat) ± 0.008(syst) (9.2) 

and infer: 

/+_= 0.490 ± 0.023. (9.3) 

Since /oo and /-| agree within experimental uncertainty, in this thesis we as- 
sume that they are equal and that the effects of any variation from (|9.ip are 
negligibly small. 

To finally estimate Aes, we compare the cut-quantity distributions of BB 
data and MC. We shift the MC distributions by small increments (the nature of 
which depend on the original distributions) to determine whether any variation 
from the original value improves the agreement between data and MC. To pro- 
vide maximum statistics, we examine all Runs simultaneously (approximately 
464 million BB events). 

If a better agreement is found, the BB efhciency of the shifted MC is cal- 
culated and the difference between this and the original value is attributed to 



Chapter 9. Systematic Uncertainty 



81 



systematic uncertainty. To measure goodness-of-fit between data and MC his- 
tograms of n bins, we use again use a statistic of the form; 

» f ATdala atMC\^ 
y^=\^ il!^ ^ ^ (9 4) 

where Nf"^"' {N^"-') is the number of data (MC) events in bin i, and for esti- 
mation purposes we make the approximations 

af*" = ^i= and af^ = — 1_. (9.5) 

The increments made to each relevant quantity are Usted below. Here, b is 
an integer, varied between and 200. 

• ETotal is shifted by an increment, up to 0.2 GeV: 

- ETotal/eCM (ETotal + 0.0015)/eCM. 

• R2 is multiplied by a constant between 0.96 and 1.04: 

- R2 ^ (0.96 + 0.0004(5) xR2. 

• To simulate the effect of tracks being lost due to tracking inefficiencies 
(when a genuine track is not correctly reconstructed as such), the tracking 
efficiency of each track (ptrack) is varied between 90% and 100%. Each 
track is given a random number between and 1, and it is kept if the 
value is less than (0.9 + 0.0005(5). Any rejected track reduces the value of 
nTracks for the event: 

- Ptrack = (0.9 + 0.0005(5). 

• The value of plMag is shifted by an increment, positive if plMag is below 
the peak of the distribution and negative if it is above the peak. This has 
the effect of narrowing the distribution: 

- plMag ^ (1.0 + 0.0005(5) xplMag if plMag < 0.97 GeV/c, 

- plMag ^ (1.0 - 0.0005(5) xplMag if plMag > 0.97 GeV/c. 

• PrimVtxdr is multiplied by a constant between 0.8 and 1.0: 

- PrimVtxdr^ (1.0 - 0.001(5) xPrimVtxdr. 

• PrimVtxdz is multiplied by a constant between 0.6 and 1.0: 

- PrimVtxdz (1.0 - 0.002^)x PrimVtxdz. 
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In each case, once the value of 5 resulting in the lowest is found (Sopt.), 
the BB MC efficiency is recalculated for the new quantity while all other cuts 
remain unchanged (ig). For some quantities like PrimVtxdz, even though the 
initial MC and BB distributions do not agree perfectly, the revised efficiency 
is almost identical to the original. This is because the cut value is so far from 
the distribution's peak that only a small number of events in the tail move 
from failing to passing the cut. The results are shown in Table IHTTl Before any 
increments are made, the original value of Sb is 

SB = 0.940456 ± 0.000033 (9.6) 

where the uncertainty stated is statistical only. 



Table 9.1: Revised efficiencies for hadronic cut quantities after shifting BB 
MC. The statistical uncertainty on each Eg value is ±0.000033 (or 
±0.003%) 



Quantity 


Sopt. 






ETotal 


87 


0.944242 


+0.403% 


R2 


110 


0.940378 


-0.008% 


Ptrack 


191 


0.940308 


-0.016% 


plMag 


8 


0.940408 


-0.005% 


PrimVtxdr 


105 


0.940644 


+0.020% 


PrimVtxdz 


59 


0.940408 


-0.005% 



The quantity with the largest change in efficiency is ETotal. The value of 
Sopt. in this case implies that adding 87 MeV to each MC event improves the fit 
to data the most. The efficiency increases by 0.403% because events up to 87 
MeV below the cut now pass selection. The original and revised distributions 
are shown in Fig. 19.31 and 19. 5i and the differences between data and MC for the 
two plots are shown in Fig. 19.41 and 19.61 

The change in R2, plMag and PrimVtxdz are negligible in comparison to 
ETotal, and the variation is comparable to the statistical uncertainty. The plots 
of data and revised MC for these variables are shown in Fig. 19.71 - 19.111 

To estimate Asb we treat each efficiency difference as an independent sys- 
tematic uncertainty and add them in quadrature. The result is: 

Aes = 0.404%. (9.7) 

9.3 Summary of Systematic Uncertainties 

We can now summarize by estimating the total B Counting systematic uncer- 
tainty for the new proposed cuts by combining the independent uncertainties 
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Total Event Energy Divided by CJUl Energy 

Figure 9.3: ETotal/eCM for BB data overlaid with BB Monte Carlo for the 
original MC ETotal cut. The solid histogram represents data (off- 
peak subtracted from on-peak) and the dotted (+) histogram is 
from MC. The histograms are normalised to have the same area 
from 0.45-1.0. 
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Figure 9.4: The difference between the data and MC histograms of Fig. 19.31 
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Figure 9.5: ETotal/eCM for BB data overlaid with BB Monte Carlo for the 
revised MC ETotal cut. The solid histogram represents data (off- 
peak subtracted from on-peak) and the dotted (+) histogram is 
from MC. The histograms are normalised to have the same area 
from 0.45-1.0. 
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Figure 9.6: The difference between the data and MC histograms of Fig. 19.51 
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Figure 9.7: R2 for BB data overlaid with BB Monte Carlo for the revised MC 
R2 cut. The solid histogram represents data (off-peak subtracted 
from on-peak) and the dotted (+) histogram is from MC. The his- 
tograms are normalised to have the same area from 0.0-0.5. 
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Figure 9.8: nTracks for BB data overlaid with BB Monte Carlo for the MC with 

revised tracking probability. The solid histogram represents data 
(off-peak subtracted from on-peak) and the dotted (-I-) histogram is 
from MC. The histograms are normalised to have the same area. 
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Figure 9.9: plMag for BB data overlaid with BB Monte Carlo for the revised 
MC plMag cut. The sohd histogram represents data (off-peak sub- 
tracted from on-peak) and the dotted (+) histogram is from MC. 
The histograms are normalised to have the same area from 0.0-5.0. 
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Figure 9.10: PrimVtxdr for BB data overlaid with BB Monte Carlo for the 

revised MC PrimVtxdr cut. The solid histogram represents data 
(off-peak subtracted from on-peak) and the dotted (-I-) histogram 
is from MC. The histograms are normalised to have the same area. 
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PrimVtxdz (cm) 

Figure 9.11: PrimVtxdz for BB data overlaid with BB Monte Carlo for the 
revised MC PrimVtxdz cut. The solid histogram represents data 
(off-peak subtracted from on-peak) and the dotted (+) histogram 
is from MC. The histograms are normalised to have the same area. 



in eb and k. Note that in Chapter [7] an optimal cutset was found by minimiz- 
ing the time- variation, measured by the expression ()7.9p . This cannot be used 
when estimating the actual systematic uncertainty in the number of B mesons, 
since by treating all MC variation as independent, it double-counts much of the 
time- variation. For the true estimate of the systematic uncertainty, we instead 
only consider the uncertainties of those quantities explicitly appearing in the B 
Counting formula, i.e. k and eb- The statistical and systematic uncertainties 
on these quantities are summarized in Table 19.21 

Combining the uncertainties (systematic and statistical) on in quadra- 
ture, we have: Eb = 0.9405 ± 0.0051. Similarly, we have the value of k = 
0.9979 ± 0.0012. The total uncertainty on the number of B mesons in a sample 
can now be found by the standard method of error propagation. From the B 
Counting formula (14.171) we can write: 



We have 



and 




With this, we are able to make an estimate on the total systematic uncer- 
tainty on the number of BB events, AA^^^. By considering the two terms of 
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Table 9.2: Summary of B Counting statistical and systematic uncertainty com- 



ponents of K^, Kx and eb- 





Uncert. 


Value 


Contrib. to AiV^ 




Statistical 


0.044% 


0.10% 




Time var. of mu-pair efF. 


0.032% 


0.07% 




Var. of 


0.047% 


0.11% 




Total 


0.07% 


0.16% 


Akx 


Statistical 


0.016% 


0.04% 




Time var. of continuum eff. 


0.096% 


0.21% 




ISR/27 contribution 


0.020% 


0.04% 




Total 


0.10% 


0.22% 


Aeb 


Statistical 


0.004% 


0.004% 




Low track-multiplicity events 


0.360% 


0.36% 




BB data/MC comparison 


0.404% 


0.40% 




Total 


0.54% 


0.54% 




Total 


0.6% 



(|9.8p separately we can determine how the uncertainties on k and Eb propagate 
into uncertainties on N'^. 

The uncertainties on and kx combined in quadrature lead to a Ak value 
of ±0.12%. This corresponds to a ±0.27% uncertainty in Ng while the Ae^ 
value of ±0.54% corresponds to an uncertainty of ±0.54% in Ng. These, added 
in quadrature via ()9.8p give the total uncertainty of 0.6% on the number of B 
mesons. 

The number of BB events counted in the BaBAR dataset is (464.8±2.8) x 10'^. 
The equivalent number given by the existing B Counting code is (465.0 ± 5.1) x 
10^. The number counted in each Run is shown in Table [9731 
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Table 9.3: Number of B mesons events with revised overall uncertainty in Runs 



1-6. 





Run 1 


Run 2 


(new cutsets) 


(22.4 ±0.1) X 10** 


(67.6 ± 0.4) X 10*^ 


Ng (existing code) 


(22.4 ±0.3) X 10'^ 


(67.4 ±0.7) X 10'^ 




Run 3 


Run 4 


Ng (new cutsets) 


(35.7 ±0.2) X 10*^ 


(110.7 ±0.7) X 10« 


N'^ (oxistiui; (:()(1(^) 


(;-!.").() ± 0.4) X 10" 


(110.5 ± 1.2) X 10" 




Run 5 


Run 6 


A^g (new cutsets) 


(146.2 ±0.9) X 10** 


(82.2 ±0.5) X 10*^ 


N]^ (existing code) 


(147.2 ± 1.6) X 10** 


(82.0 ± 0.9) X 10** 




Total 




TV^j (new cutsets) 


(464.8 ±2.8) X 10« 




Ng (existing code) 


(465.0 ±5.1) X lO'* 
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Chapter 10 



Summary 



We have presented a new proposed set of cuts designed to improve the precision 
of B Counting at the BaBAR Experiment. The overall uncertainty in the number 
of B mesons counted has been reduced from 1.1% to approximately 0.6%. 

BGFMultiHadron and nTracks are the largest contributors to uncertainty 
due to inaccuracies in simulating low track-multiplicity events 0.36%) and 
from the differences between the ETotal distributions of BB data and MC (~ 
0.40%) . 

Algorithms implementing these new B Counting selection criteria are being 
run in the Release 24 reprocessing of the BaBAR datasct. This ensures analysts 
running on the most recent version of the data will be able to count B mesons 
with this improved method. To allow cross-checking, the old tagbits isBCMuMu 
and isBCMultiHadron are also recalculated during reprocessing. 
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Uncertainties of a 
Multi- Variable Function 



Let f{xi,---,Xm) be an infinitely differentiable function of m variables {xi), 
each with an uncertainty Axi. The Taylor series of / expanded around the 
point (ai, - • • ,am) is: 

T(x.---x ) - V V ^""^ f{au---,am) 

•(a;i-air---(a;„-a„r-. (A.l) 

Assume each Xi comes from an independent statistical distribution with 
mean Xi and standard deviation Ax^ . If the standard deviations are sufficiently 
small, we can approximate / by evaluating the series to first order around 
(xi, • • • jXm). The uncertainty of the function can then be estimated by the 
equation: 

[Af{x^,---,xm)f =Y,^(^^Ax}j , (A.2) 

where the differential is evaluated at (xi, ■ ■ ■ ,Xm)- The zeroth order terms are 
constant and do not affect the uncertainty. 

The expression for /t^^ has the same form as 

/(xi,X2) = ^, (A.3) 

SO in this case: 

(") 

Similarly, Kcont. has the form 



and 



2 



g{xi, X2, xs, X4) = "^^^ (A.5) 

X3 + Xi 



Axi \ ( A.T2 \ f -{Xl + X2)AX3, 



xs + Xi J \ is + X4 y V XiY 

Xi + ^2)^ 

(f 3 + X4,y 
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Appendix B 

T{3S) Counting 



Between 22 December 2007 and 29 February 2008 (Run 7), the BaBAR detector 
recorded events provided by PEP-II at the T(35) resonance. Data was provided 
at the resonance peak with a CM energy of 10.35 GeV (on-peak data) and also 
at approximately 30 MeV below this (off-peak data). 

The author and Chris Hearty performed an analysis of T{3S) Counting on 
this dataset. The methods used were closely based on the work presented in 
this thesis and [11] , and are described in detail in a BaBAR Analysis Document 
[28]. 

During Run 7, 121.9 million Y{3S) mesons were produced. These were 
counted to within a systematic uncertainty of 1.0%. Some results utilizing this 
work have been published, including a significant observation of the bottomo- 
nium ground state r]i,{lS) through the decay T(3S') — > jr]b{lS) [29] . 
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